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Abstract. We study questions related to critical points of the Green's func- 
tion of a bounded multiply connected domain in the complex plane. The 
motion of critical points, their limiting positions as the pole approaches the 
boundary and the differential geometry of the level lines of the Green's func- 
tion are main themes in the paper. A unifying role is played by various affinc 
and projective connections and corresponding Mobius invariant differential op- 
erators. In the doubly connected case the three Eisenstein series E2, E4, Eg 
are used. A specific result is that a doubly connected domain is the disjoint 
union of the set of critical points of the Green's function, the set of zeros of 
the Bergman kernel and the separating boundary limit positions for these. At 
the end we consider the projective properties of the prepotential associated to 
a second order differential operator depending canonically on the domain. 
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1. Introduction 

The results of this paper have their origin in attempts to understand trajectories 
of critical points in a planar multiply connected domain. We start with study- 
ing, for a given bounded multiply connected domain f2 in the complex plane, the 
motion of critical points of the ordinary Green's function G(z, £) of as the pole 
£ moves around. If Q has connectivity g + 1 and the critical points are denoted 
Zi(C)j ■ ■ ■ , %(C)j we show first of all that the Zj(() stay within a compact subset K 
of Q. (this result has also recently been obtain by A. Solynin [56]) and secondly that 
the limiting positions of the Zj(Q as C approaches the boundary coincide with the 
corresponding limiting positions of the zeros of the Bergman kernel. In the doubly 
connected case (g = 1) it even turns out that the domain is the disjoint union of 
the set of critical points of the Green's function and the set of zeros of the Bergman 
kernel, plus the common boundary of these two sets. 

The method developed to prove the existence of the compact set K is remarkably 
related to (a new type of) Martin compactification and uses many properties of the 
Bergman, Schiffer and Poisson kernels. All these functionals depend in one sense 
or another on critical points of the Dirichlet Green's function, and they moreover 
extend in a natural way to sections of suitable bundles of the Schottky double of the 
domain, which is a compact Riemann surface of genus g. An important role here is 
played by the Schwarz function, which can be interpreted as being the coordinate 
transition function between the front side and the back side of the Schottky double. 

It is natural in this context to study the relationship between the hyperbolic 
metric and level lines of the Green's function. We show, for example, that the 
average, taken on a level line of the Green's function, of the Green's potential 
of any compactly supported measure /i is proportional to the total mass of /i. We 
also study level lines of the Green's function, and other harmonic functions, from 
the point of view of hamiltonian mechanics and differential geometry. One main 
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observation then is that such level lines are trajectories, as well as Hamilton- Jacobi 
geodesies, for systems with the squared modulus of the gradient of the harmonic 
function as hamiltonian. 

Various other topics touched on are estimates of the Taylor coefficients of the reg- 
ular part of the complex Green's function in terms of the distance to the boundary, 
and relationships between these coefficients and the Poincare metric. 

A major part of the paper is devoted to a detailed study of the doubly connected 
case, modeled by the annulus. In this situation all calculations can be done explic- 
itly because we have at our disposal the theory of modular forms (theta functions) 
and elliptic functions. For example the dichotomy result mentioned in the begin- 
ning is proved in this section. The Ramanujan formula for the derivatives of the 
basic Eisenstein series £2, E4, Eq are fundamental. 

Some of the mentioned Taylor coefficients and Eisenstein series transform under 
coordinate changes as different kinds of connections. In the second half of the 
paper we study affine and projective connections in quite some depth. In general, 
projective connections play a unifying role in the paper, in fact almost all of our 
work turns around projective connections, and to some extent affine connections. 
One example is that they generate Mobius invariant differential operators, called 
Bol operators and denoted A m , for which boundary integral formulas, of Stokes 
type but for higher order derivatives, can be proved. These are useful for studying 
weighted Bergman spaces in general multiply connected domains. There is one 
such Bergman space for each half-integer, and the elements of the space should be 
thought of as differentials of this order (actually integrals if the order is negative). 
The Bol operator then is an isometry A m : A m (Q) — > B m (Cl) from a Bergman type 
space A m (Q) of differentials of order i^ 21 to a (weighted) Bergman space B m (Q) 
consisting of differentials of order -^-j 2 . Also reproducing kernels for these, and 
some other, spaces are studied. 

Along with connections, several questions related to curvature of, for example, 
level lines are discussed. The Study formula finds here its natural meaning. A 
further aspect is that the operators A m all turn out to be symmetric powers of a 
single one, namely A 2 . 

In the final section of the paper we try to clarify in our setting the meaning of a 
certain prepotential for second order differential operators which has arisen in some 
recent physics papers. 

2. Generalities on function theory of finitely connected plane 

DOMAINS 

2.1. The Green's function and the Schottky double. Let f2 C C be a bounded 
domain of finite connectivity, each boundary component being nondegenerate (i.e., 
consisting of more than one point). The oriented boundary (having f2 on its left 
hand side) is denoted dil = T = r + ri + - • - + r g , where r is the outer component. 
We shall in most of the paper discuss only conformally invariant questions, and 
then we may assume that each boundary component Tj is a smooth analytic curve. 
Alternatively, one may think of Q as a plane bordered Riemann surface. 

Much of the functions theory on O is conveniently described in terms of the 
Schottky double Cl of il. This is the compact Riemann surface of genus g obtained, 
when Q has analytic boundary, by welding f2 along 9f2, with an identical copy 0. 
Thus, as a point set = ft U dfl U ti. The "backside" f2 is provided with the 
opposite conformal structure. This means that if z G denotes the point opposite 
to z G Q then the map z 1— ► z is a holomorphic coordinate on Q. The construction 
of the Schottky double generalizes to any bordered Riemann surface and the result 
is always a symmetric Riemann surface, i.e., a Riemann surface provided with an 
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antiholomorphic involution. In the Schottky double case this is the map J : Cl — > Q 
which exchanges z and z and which keeps dQ. pointwise fixed. 

The Schottky double of a plane domain £1 has a holomorphic atlas consisting 
of only two charts: the corresponding coordinate functions are the identity map 
4>i : z i— ► z on Q and the map <p2 '■ z i— > z on f2. When dQ is analytic, as is 
henceforth assumed, both these maps extend analytically across dfl in Q, hence 
their domains of definitions overlap and the union covers all f2. Let 

(2.1) S = 2 o 

be the coordinate transition function. It is analytic and defined in a neighbourhood 
of dfl in C, and on d£l it satisfies 

(2.2) S(z) = z (ze dfl). 

Thus it is the Schwarz function [TU], [H], [53] of d£l. 
Differentiating (|2.2|) gives 

(2.3) dz = S'(z)dz along dfl. 

With s an arc-length parameter along dil such that f2 lies to the left as s increases, 

(2.4) TM-* 

is the oriented unit tangent vector on dil. By (12. 3p and since |T(z)| = 1, 

It follows that 1/T(z) extends analytically to a neighbourhood of dil and that it 
gives a selection of a square- root of S'(z). This is an important observation because 
it means that on the Schottky double of any plane domain there is a canonical 
choice of square-root of the canonical bundle, i.e., there is canonical meaning of 
the concept of a differential of order one-half, and hence of a differential of any 
half-integer order. 

A function / on Q is most conveniently described as a pair of functions /i, f% 
on f2, continuously extendable to dil, such that 



f 1 (z) = f 2 (z) (zedty. 

The formal relations to / in terms of the coordinates functions <j)\ and <f>2 above are 

\h = co f ocj)- 1 oc, 

where c denotes complex conjugation. It follows, for example, that / is meromorphic 
if and only if fi and f% are meromorphic. With the same kind of identifications, a 
differential of order m on il is represented by a pair of functions f±, f% on f2 such 
that 

fi{z)dz m = h{z)dz m along dQ. 
This is to be interpreted as fi(z)T(z) m = f 2 (z)T(z) m , or 



(2.5) h{z) = h{z)T{zf™ (zeflfi). 

Clearly this makes unambiguous sense for any m € iZ. 

The Green's function G(z, Q of Q is, as a function of z for fixed £ S fi, defined 
by the properties 

G(z, C) = — log \z — C| + harmonic zell, 
G(z, () = 0, zedn. 
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It is symmetric in z and £, G(z, £) = G((, z), and it extends to the Schottky double 
as an "odd" function in each variable, for example, G(J(z),() = —G(z,Q. 

The above extension of the Green's function makes it a special case of a fun- 
damental potential which exists on any compact Riemann surface, and which is a 
suitable starting point for discussing the classical function and differentials: for any 
three distinct points a, 6, w on a compact Riemann surface M there exists a unique 
function of the form 

(2.6) V(z) = V(z, w; a,b) = — log \z — a\ + log \z — b\ + harmonic, 
normalized by V(w, w; a, b) = 0. It has the symmetries 

(2.7) V(z, w; a, b) — V(a, 6; z, w) = —V(z, w; b, a) 
and the transitivity property 

(2.8) V{z, w; a, b) + V(z, w; b, c) = V(z, w; a, c). 

See below for explanations, and also [52], Ch. 4. If The Riemann surface is sym- 
metric, with involution J, then 

(2.9) V(z, w; a, b) = V(J(z), J(w); J (a), J(b)). 

Example 2.1. In the case of the Riemann sphere, M — P, we have 

[z — a)(w — b) I 



V(z, W] a,b) = — log \(z : w : a : b)\ = — log 
(z : w : a : b) denoting the classical cross-ratio. 



{ z -b){w- a y 



With M = fl, the Green function is given in terms of V by 

(2.10) G(z, = ±V(z, J(z);(, J(C)) = V(z, w; (, J(Q), 

where w is an arbitrary point on dtt and where the second equality follows from 
(I2T7D . (23J, (EH)- Cf. also HE], p. 125f. 

The existence of V(z, w; a, b) in general is immediate from classical potential 
theory, see e.g. [S2], [H] or, more generally, from Hodge theory. In fact, V solves 
the Poisson equation — d*dV — 27r(<5 a — (5b) on M, where the star is the Hodge 
star and the Dirac measures in the right member are to be considered as 2-form 
currents; the solution exists because J M 2n(S a — 5b) = and it is unique up to an 
additive constant. 

From V(z, w; a, b) much of the classical function theory on M can be built up. 
For example, 

(2.11) v a _ b (z) - -2 dV{z '™ ;a ' b) dz = -dV{z) - i*dV(z) 

oz 

dz dz . 

analytic 



z — a z — b 

is the unique abelian differential of the third kind with poles of residues ±1 at 
z = a,b and having purely imaginary periods. The subscript in u a _b should be 
thought of as a divisor, and the definition extends to any divisor of degree zero. 
Conversely, we retrieve V from v a ^ b by 

(2.12) V(z,w;a,b) = -Re f v a _ b , 

J w 

from which the symmetries (|2.7p and transitivity property (|2.8p of V follow, using 
Riemann's bilinear relations for one of the symmetries. 

From V(z,w;a,b) one can also construct the 2g-dimensional space of harmonic 
differentials on M by considering the conjugate periods. By (|2.12p . the harmonic 
differentials dV = —Kev a -b and *dV = — lmv a -b do not depend on w. The first 
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one is exact, while the second has certain periods, which depend on a, b: if 7 is a 
closed curve on M then the function 

<t ) y( a , b ) = 7T- / *dV(-,w;a,b) 

Z7T ' 



7 



is, away from 7, harmonic in a and b and makes a unit jump as a (or b) crosses 7. 
It follows that, keeping b fixed, d(f) 1 (- 1 b) extends to a harmonic differential on M 
with periods given by 

/ d(j) 7 {-,b) = a x 7, 

J a 

where a x 7 denotes the intersection number of cr and 7. Choosing 7 among the 
curves in a canonical homology basis gives the standard harmonic differentials. 

As a further aspect, the right member in (|2.12|) can be written as a Dirichlet 
integral: 

(2.13) V(z,w;a,b) = ~ [ dV(-, c; z, w) A* dV(-, c; a, b) 

2tt Jm 



m 

= — Im / Uz-w A -D a _ & , 

where c € M is an arbitrary point. Thus V reproduces itself in a certain sense, 
and the equation also expresses that V(z, w; a, b), besides being the potential of the 
charge distribution 5 a — 6b when considered as a function of merely z, also is the 
mutual energy of the two charge distributions S a — 6b and 6 Z — 6 W . The first equality 
in (|2.13p follows by partial integration, and the second from the definition of v a -b- 
For the Green's function the corresponding equation is 



G(z,Q = ^- [ dG(-,z) A* dG(-,0- 
27r Jn 



We recall that v a -b has 2g zeros (since it has 2 poles). With M = & and 
a = C, b — J(C), half of the zeros are on f2, and by (|2.10p . (|2.1ip these are exactly 
the critical points of G(z, £) (the points where the gradient vanishes). Hence the 
Green's function has exactly g critical points on Q. 

2.2. The Schottky-Klein prime function. The harmonic theory on a compact 
Riemann surface is simple and intuitive, as indicated above. For the holomorphic 
theory one usually prefers to work with abelian differentials of the third kind nor- 
malized, not as v a —b, but so that, in terms of a canonical homology basis, half of 
the periods vanish. Following [SS] and [3S] we shall, in the case of a Schottky double 
of a planar domain fi, choose a canonical homology basis {ai, . . . , a g , j3\, . . . , /3 g } 
such that ctj goes from To to Tj on f2 and back to To along the same track on the 
back-side f2, and such that — Tj (j = 1, . . . ,g). Thus the basis is symmetric 
with respect to the involution J, more precisely J (ctj) — —ctj, J((3j) — /%. 

We denote by aj a -b the abelian differential of the third kind with the same 
singularities as v a -b but normalized so that the a^-periods vanish: 

(2.14) f uj a _ b = 0, i = l,...,g. 

This makes sense only for a,b (£ ctj, and with the ctj being fixed curves. Hence 
LJ a -b is less canonical than v a —b but it has the advantage of depending analytically 
on a, 6, whereas for v a -b the dependence is only harmonic in general. 

Remark 2.1. As to notation, the differential we denote by w a -6 is in Fay [18] . 
Schiffer-Spencer [55], Farkas-Kra [T7] denoted, respectively, u> a -b, —duj a b, T a b- For 
our v a -b the corresponding list is: fl a -b, —d£l a b) ^ab- Also, Fay [TH], and Hejhal [55] 
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use a homology basis with switched roles between the aj and (3j curves, but, as no- 
ticed by A. Yamada [55], the present choice has certain advantages (see Lemma HOI 
below). 

The integral f* ui a -b, with unspecified path of integration, is locally holomorphic 
in all variables, but multivalued. To cope with the 2iri indeterminacy coming from 
the poles one may form the exponential. In the case of the Riemann sphere this 
simply gives the cross-ratio between z, w, a, b: 

exp / uj a -b = (z : w : a : b). 

J w 

For Riemann surfaces of genus g > 0, the exponential remains multivalued, with 
multiplicative periods (cf. [17]), Ch.III.9). However, one can still write it as a kind 
of cross-ratio if one is willing to accept further multivaluedness: 

T E(z,a)E(w,b) 

(2-15) exp / u a _ b = -. 

Jw E(z,b)E(w,a) 

Here E(z,£) is the Schottky-Klein prime function (prime form), which should be 
regarded as a differential of order — ^ in each of the variables, but as such still has 
multiplicative periods. Near the diagonal it behaves like z — £. The exact definition 
of E(z, is usually given in terms of theta functions on the Jacobi variety, see 
[18] , [28] , [42] and below. For the Schottky doubles of a plane region there is also 
a representation in terms of the Schottky uniformization of the domain, see for 
example [Hi- 
lt should be remarked that, in the Schottky double case, the kind of differential 
of order — \ referred to for E(z, £) has no representation in the form (|2.5p (with 
rn = —1). The bundle of half-order differentials defined by (|2.5[) is "even" and 
does not allow any global holomorphic section, whereas the one needed for E(z, Q 
should be "odd" , which does allow for a holomorphic section. The definition of 
such a bundle requires a finer atlas than the one consisting of only O and O for its 
representation. Of the 2 2g bundles of half-order differentials, 2 s ~ 1 (2 g — 1) are odd 
and 2S- 1 (2« + 1) are even, see [28]. 

Now, quite remarkably, in case the Riemann surface is the Schottky double of 
a plane domain the two abelian differentials Lu a -b and v a -b coincide when a and b 
are symmetrically opposite points: 

Lemma 2.1. For M = Cl and with canonical homology basis chosen as above, 

Remark 2.2. The lemma does not hold if the homology basis is chosen in a different 
way, for example with switched roles between the ctj and (3j curves. Note also that 
even though u) a —b makes a jump by 27ri as a or b crosses aj, there is no such jump 
for ujq_j(q because C and J(C) cross aj simultaneously and the two contributions 
cancel. 

Proof. One simply has to notice that v^_j(q satisfies the normalization (|2.14p of 
W C— J(0' E x P resse d in terms of V, since dV already is exact, the statement to be 
proven becomes, for any w £ dfl, 

*dV(;w;C,J(0)=0, (j = l,...,g). 

That these periods vanish follows from the symmetry properties (|2.13[) . (|2.9p of V 
with respect to J. □ 
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For a general Riemann surface, the differential oj a -b can be recovered from the 
prime form as a logarithmic derivative of (|2.15[) : 

(2.16) c Q _ b (z)=dlog||^|. 

From this one sees clearly that E(z, £) depends on the homology basis, since Lo a -b 
does. The Green's function of a planar domain Q is most directly related to w a -b 
on ft as in (|2.10p . (|2.1ip . but in view of Lemma |2~T1 one can equally well use w a -&- 
Therefore, (|2.16[) gives the following expression of the Green's function of SI in terms 
of the prime form on M = Q, 

G ^ c ) = - l0g l^7§))l ^ e0 >- 

See 59 for further details, and also 9 , 35 for other aspects. 

We introduce next the harmonic measures Uj, j — 1, . . . ,g, i.e., the harmonic 
functions in f2 defined by having the boundary values Uj — 5k j on Tfc. It is easy 
to see that their differentials duj extend to the Schottky double as everywhere 
harmonic differentials with 



duj = 25k j, / *duj = 0. 

*ot k Ja k 

Here the second equation is a consequence of the symmetry of the extended dif- 
ferential under J. Thus, the dlij — \{duj +\*duj) (j = l,...,g) constitute the 
canonical basis of abelian differentials of the first kind (everywhere holomorphic 
differentials), the period matrix with respect to the a& curves being the identity 
matrix. For the /3fe-periods we have 



duj — 0, / *duj — Uk *duj = / duk A *duj, 
p k J/3 k Jan Jn 

where the latter make up a positive definite matrix. On setting 

Tkj = ~ / duj +i*dv,j 
1 Jh 

we thus have ReT^- = and that the matrix (Imr^ ) is positive definite. With 
Uj{z) = J z \{duj +i*duj) = f z ^-dz denoting the corresponding abelian integrals 
(multivalued), the map 

z^U(z) = {Ui{z),...U g {z)) 

defines, up to a shift, the Abel map into the Jacobi variety C s /(Z s +rZ g ), r = (r^-j). 

At this point we can make the definition of the Schottky-Klein prime function 
slightly more precise. The first order theta function with (half-integer) character- 



istics 
(2.17) 

i) 



where 5, e £ \V are row vectors, is defined by 



'5 


(w) = & 


s 






c 




c 



(w; t) = exp [m(m + 5)r(m + + 2tti(w + e)(m + Sf] 

for w — (wi, . . . ,w s ) G C s . The superscript t denotes transposition of a matrix. 
To define the prime form the characteristics should first of all be chosen to be 
odd, i.e., so that 4Se t is an odd number (e.g., 5 — (|, . . . , 0), e = (^, . . . , 0)). 

whenever such a choice of 5, e has been made, and 



For simplicity, set i?* = d 
introduce 
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The S and e should in addition be chosen to be non-singular, i.e., so that not all of 
the constants Cj vanish. This can always be done, see [TB]. Then the Schottky-Klein 
prime form is, when considered as defined on Cl, 

E , c) = fl.(M(z)-M(0) 

See 35 (Appendix) for further details using the present point of view. 

3. Critical points of the Green's function 

3.1. The Bergman and Poisson kernels. The Bergman and Schiffer kernels [3], 
[52], K(z,Q an d L(z,Q respectively, will be discussed in some detail later on. For 
the moment we just recall their representations in terms of the Green's function 
and the potential V(z, w; a, b) in (|2.6p : 

(3-1) K(z,Q = --~(z,Q= -~(z,J(z);QJ(Q), 

7T OZOQ 7T ozoQ 



By the symmetry of G{z, Q (or V) we have L(z, Q — L(Q z) and K{z, Q — K(Q z). 
Since G(z, Q = for z G dSl, ( G Q, 



(3.3) L(z, Qdz + K (z, Qdz = 

along dfl (with respect to z). This means that, for any fixed £, the pair K(z, Qdz, 
L(z, Qdz combines into a meromorphic differential on il. It has a pole of order two 
residing in L(z, Q: 

(3-4) L(z,Q = - r ^ 7 ^-i(z,Q, 

where £(z, Q (the "^-kernel") is regular in both variables. Consequently, K(z, Qdz, 
L(z, Qdz have altogether 2g zeros. For £ G c?S7 the zeros are, by (|3.3[) with switched 
roles between z and £, equally shared: g zeros for K(z, Qdz and L(z, Qdz at the 
same points, and no zeros on <9f2 as we shall see in the next section. A further 
consequence of (I3.3[) is that when <9f2 is analytic, the kernels K(z,Q, L(z,Q have 
analytic extensions across dil in C. 

Using the symmetries of K(z, Q and L(z,Q, and (|3.3|) one finds that the double 
differentials K(z, QdzdQ L(z, Qdzd( are "real" on the boundary: 

L(z,Qdzd( e M, K(z, QdzdC, G M (^,(6 90). 

The precise meaning of for example the first equation is that L(z, QT(z)T(Q G R, 
where T(z) is the tangent vector at z G 9£1 (see (12.41) ). 

The Poisson kernel of fl is the normal derivative of the Green function when one 
of the variables is on the boundary: 

P(z,Q = -——(z,Q, zedn,(en. 

2ir on z 

Here denotes the outward normal derivative at z G With ds = ds z denoting 
arc-length with respect to z, the definition may also be written in differential form 
in several ways, for example (with d = d z ) 

P(z,Qds z = -J-*dG(z,C) = -~Im^-(z,Qdz 

Z7T 7T OZ 

1 dV 1 
= —^ l ^^(z,J(z);QJ{Q)dz = -— Im« ( _j( ( )(4 
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3.2. Critical points. In this subsection and the next we show that, as the pole 
moves around, the set of critical points of the Green's function of fi stay within 
a compact subset of fi, and that the limiting positions, as the pole tends to <9fi, 
coincide with the zero set of the Bergman kernel in this boundary limit. The latter 
result was the main motivation for the present work. Related results have been 
obtained by S. Bell [3] (section 30) and A. Solynin [55] , 

For fixed z € <9f2, the map £ i— > P(z, C) is harmonic and strictly positive in 
51 and vanishes on dfl \ {z}. It extends, when dfl is analytic, to a harmonic 
function in a neighborhood (in C) of dfi except for a pole at £ = z. Since P(z, •) 
attains its minimum in Sl\{z} on <9fl, the Hopf maximum principle 47 shows that 

|^-(z, C) < for all C e dQ \ {z}. Thus 

(3.5) J^L( Z , C ) >0 , z,(effl, z^C- 

on z on^ 

The fact that this double normal derivative is nonnegative can be more directly 
understood by interpreting the left member of (|3 . 5|) as a double difference quotient: 
with differences Az, AC pointing in the normal direction into il from z,( e 9J1, 
respectively, we have 

d 2 G(zX) _ 
dn z driQ 

G(z + Az, C + AC) - G(z + Az, () - G{z, ( + AC) + G(z, 

— J- lm n. — rmr, • 

|A*|,|AC|->o |Az||AC| 

This is obviously nonnegative because all terms in the numerator vanish except the 

first: G(z + Az,C + AC) > 0. 

Complementary to (|3.5p we have 

d^G . ,_. d^G , _. O^G . , , 

It follows that, for £ E 90, 2 7^ C 

O^G , . _ O^G , , „ d^G , ,_•./„ 

^c ( ^ ' ^C ( ^ ' ^ ( ^ ' 

because these derivatives contain derivations in the normal direction in each of the 
variables. In view of the singularity of type l/(z — C) 2 at z = C and of the continuity 

of and = as mappings from a neighborhood of 951 x <9£! c C x C to the 

ozoQ ozd( 

Riemann sphere it is clear that the above inequalities persist to hold for z = £, and 
also that the quantities are bounded away from zero: 

in a full neighborhood of dil x <9il and for some constant c > 0. Note that (|3.6p says 
that -fC(z,C), L{z, C) are bounded away from zero in a neighborhood of 911 x <9£1 
when c?fl is analytic. With z in a neighborhood of 9J1 and C € 9f2 we also infer 

(3.7) |^^|> c>0 . 

1 az oriQ ' 

dG 

Since —— (z, O = for C G c?fl we obtain, on integrating (13. 7D with respect to C 
oz 

in the normal direction, 

BC 

| — (z,C)| > cdist(Cdfi) 
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for z, £ in a neighborhood of dQ and for some c > 0. In particular there exists a 
compact set K C such that 

(3.8) ^-teO^o 

oz 

for all z, C € £1 \ K. Now we are ready to conclude the following, also obtained 
(using slightly different arguments) by A. Solynin [56 . 

Theorem 3.1. Let Q C C be a bounded finitely connected domain such that each 
component of C \ f2 consists of a least two points. Then there exists a compact set 
K C f2 such that 

dC 

(3.9) a^ z ' C ) *°> * e °\^> Cefi, 

(3.10) S| (z ' C) *>Cen\#, 

(3.H) ^0, z,(€J2\if. 

Proof. The statements are conformally invariant, so it is enough to prove them when 
O has smooth analytic boundary. By (|3.6[) . (|3.7[) . (|3.8p the desired inequalities are 
valid for z,( e Q \ K and it remains only to prove that, in the first inequality, we 
can allow all £ £ f2 by possibly enlarging the compact set K. If this were not true, 
then there would exist a sequence {(z n , Cn)} with 

QQ 

QQ 

According to l|3.8[) . ( e U. But then also (z,£) = 0, which however cannot 

QQ 

be true because (z, C) is a nonzero constant factor times the Poisson kernel 

dz 

P(z, C), which is strictly positive for £ £ Q, z S dfi. This contradiction finishes the 
proof. □ 

Remark 3.1. With slight modifications the proof also works in several real variables: 
if £1 C M. n has real analytic boundary then there exists a compact set K C ri such 
that Va;G(a;, ?/) ^ for all x £ fi \ if , y £ ft. The crucial observation is that (|3.5[) 

n d 2 G(x y) 

persists to hold, so that — — — - — ^rjj ^ for x, y £ dQ whenever £, r\ are 
nontangential vectors at x and y respectively. 

Example 3.1. The boundary behavior and the nature of the pole are illustrated by 
the case of the unit disk: = D = {z £ C : |z| < 1}. In this case Cl = P = 
CU {oo} with involution J : z i— > 1/z, and the Schwarz function is S(z) — l/z. The 
fundamental potential was given in Example 12.11 and the Green's function is 

■c, 



It follows that 



G(*,C) = -log| 

l-z( 

dG(z,0 _ 1-|C| 2 



dz 2(z-C)(l-zC) : 
d 2 G 11 d 2 G 1 



9z9C 2(z-C) 2 ' 9z9( 2(l-zC) 2 

i-ICI 2 



^,0 = 



CI 2 
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As for the estimates (|3.6[) . p.7[) we note that \z — C| < 2, |1 — zC\ < 2 so that, for 

, d 2 G , 1 , d 2 G , 1 , d 2 G , 1 
\dzdC ~ 8' Uz^c' ~ 8' '<9zdn c ' ~ 8' 

3.3. Critical points in the boundary limit. When C £ c?S7, — — (z,C) = for 

dG 

all z G fi, but by "blow-up" one can still speak of a nontrivial limit of — — (z,C) 

oz 

d 2 G 

as C — > d£l. In a certain sense 7— - — (2, C), and hence the Bergman and Schiffer 

ozon^ 

kernels, represents this limit, but there is also a representation in terms of a Martin 
type construction, which we shall discuss in the next section. Throughout this 
subsection we assume that <9f2 is analytic. 

Fix a point a 6 f2 such that d z G(a, () ^ for all C G fi and if (a, C) 7^ for all 
C S dfi. This is possible (with a close enough to the boundary) by Theorem 13.11 
and what precedes it, viz. (|3.6p . Then define 

(3.12) FM= WM zXe n. 

d z G(a,Q 

As a function of z, F(z, £) is meromorphic in SI with a pole at £. It is normalized so 
that F(a, C) = 1, which prevents it from degeneration as £ — > cT2. In Example 13.11 
for the unit disk, we see this from 

( a -Q(l- g Q _ (0-0(0-1/0 

1 ' u (2-c)(i-<) (z-0(*-i/CT 

As a comparison, 

if(z,Q = (l-<) 2 = (a-l/C) 2 
tf(a,C) (1-zC) 2 (^1/C) 2 ' 

Since £ = l/C on the boundary, it follows that F(z, £) and — ; ' ^ coincide there. 

AT(a,C) 

The following theorem shows that the above is what happens in general. 

Theorem 3.2. Assume dfl is analytic and choose a £ £1 as above. Then the 
function F{z,Q, originally defined in fi x ft, extends continuously to ft x ft, and 

K(z,() 

onSlx dfl it agrees with — ' . Specifically, if (z„, £ n ) 6 x O, (z, G O x 30 

K(a,C) 

and (z n , Cn) - ► (z, C) flSK-t oo, £/ien 

lim F(z n ,C n ) = §f^-. 

Proof. Let (z„,£„) — > (z,0 as in the statement and let r;„ £ dft be closest points 
on the boundary to the 

|Cn ~Vn\= d(( n ) = d(C n ,dQ). 
With T(r] n ) the tangent vector at r) n £ dft, we have 

■^-(Zn,?7nJ = 0, 

<9 2 G , , „, s <9 2 G 



(z n ,r] n )T(rj n ) + T-7^(z n ,rj n )T(r] n ) = 0. 



ftrflC' " ty 9z9C 
Therefore, Taylor expansion with respect to £ of ^-(z„, £) at £ = ^„ gives 

<9G <9 2 G <9 2 G 

T-(Zn,Cn) = ^-^:(^,??n)(Cn ~ »7n) + ^"^r (z n , ??n ) (C« ~ Vn) + C(ICn ~ ??n| 2 ) 

oz ozoQ ozoQ 
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= 2i S| (Zn ' Vn) T(rin)Im + °( d ^ 

= 2ip?-(z n ,7 ln )T(r ln )d(( n ) + 0(d(( n f). 
ozoQ 

Note that Im = d(( n ). The remainder term 0(d(( n ) 2 ) in principle depends 

on z n , but it is uniformly small as z n — ► z 6 ft. 
We now obtain 

F d z G(z,C») = 2i0(z„,77»)rfan)^(Cn)+g(rf(Cn) 2 ) 

^(z»,?7») + 0(d(Cn)) 

*T(a,rM) + 0«Cn)) ' 

Since if (a, £) is bounded away from zero for ( close to dft the statements of the 
theorem follow. □ 



As a consequence we have 

Corollary 3.3. The limit set of the set of critical points of G(z, £) as C — > <9ft is 
exactly the set of zeros of K(z, C) for C G <9ft. More precisely, for z £ ft, C G <9ft; 
i/ie following statements are equivalent. 

(i) ff(*,c) = o. 

(ii) L(2,C) = 0. 

(hi) There exist (z n , Cn) G ft x ft, (n = 1, 2, ■ • • ) wii/i d z G(z n , Cn) = suc/i i/ia£ 

(^n, Cn) -* (z, C) as oo. 
(iv) For eac/i sequence {Cn}n>i C ft swc/i iftai Cn — > C> i/iere exist z„ € ft wzt/i 

z„ — > z sztc/i i/iai d z G(z n , Cn) = 0. 

Proof. By p.3[) (with switched roles of z and (), (i) and (ii) are equivalent, and 
clearly (iv) implies (iii). 

If d z G(z n ,( n ) = 0, i.e., F(z n ,(n) = 0, then it is immediate from Theorem 13.21 
that K{z,C t ) = 0. Thus (iv) implies (i). Conversely, assume for instance that 
K(z,Q = 0, ^K(z,Q ^ and let Cn € ft, Cn -> C- The functions F(-,C„) and 
if(-,C) have equally many poles and zeros (namely 2 poles and 2g zeros seen in 
the Schottky double), and since F(-, Cn) tends to a constant times K (•, C), F(-, Cn) 
must have a zero z„ near the zero z of K (•, C)- In fact, let 7 = <9B(z, e) with e > 
small enough. Then, as n — > 00, 

^/^^(•=Cn)-^/^l0g^- ) C) = l, 

the last integral being de number of zeros of K(-,Q inside 7. Hence F(-,(n) has 
exactly one zero z n inside 7. Letting e tend to as n tends to 00, we conclude that 
there are zeros z n with z n — > z. Thus (i) implies (iv), and the proof is complete. □ 

Remark 3.2. The following example illustrates the fact that critical points of Green's 
function are not necessarily simple. Let u = e 2 ^ 1 ^ 3 , Dj — D(cj j , p), j = 1, 2, 3, and 

< p < The Green's function G(z) of ft = C \ {D 1 U D 2 U As) with pole 
at infinity has two critical points (g = 2). Since ft is invariant under z — > ujz, 
the two critical points must be a double point located at the origin. By means of 
a Mobius transformation, one obtains an example of a 2-holed disk which has a 
Green's function with a multiple critical point. 
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4. The Martin and gradient boundaries 

4.1. The Martin compactification. In this section we shall give the preceding 
considerations their right meaning. We adapt the construction of Martin compact- 
ification as presented for instance in I.S. Gal [2T]. Let us first recall a general 
theorem of Constantinescu-Cornea (see [TJ, p. 97). 

Theorem 4.1. Let be a non-compact locally compact Hausdorff space and let <& 
be a family of continuous functions Q — > [— oo, +oo]. Then there exists a compact 
topological space £l M , unique up to homeomorphisms, such that 

(i) is an open and dense subset of £l M . 

(ii) Every f G $ can be extended to a continuous function f M on Q M . 

(iii) The functions f separate points on d M Vl := f2 M \ il. 

For example, when fl is a multiply connected domain in the plane one may, 

for a fixed a € O, consider the family <f> of functions C >— > M(z, £) = ^ ' - - , 

G(a,C) 

parametrized by z € f2 and with the convention that M(a, a) = 1. Each function 
z i ► Af (z, C) is continuous, even at a. The space f2 M obtained by the theorem 
of Constantinescu-Cornea is, up to a homeomorphism, independent of a. It is the 
Martin compactification of f2, and the Martin boundary is d M VL = VL M \ Q, 

As discovered by M. Brelot and I. S. Gal (see [21]), the Martin boundary can 
also be introduced via uniform structures. For generalities on such, see |32j . Let 
X be a nonempty set, let (Y, V) be a uniform space and let $ be a family of 
functions <f> : X — > Y. Then there is a weakest uniform structure on X making the 
functions in $ uniformly continuous. It is the uniformity U generated by all the 
0- x (V) = {{x,y) e X x X : [<j>{x) , (j>{y)) € V}, for V £ V, <f> £ $, as a subbase. 

If V is precompact (totally bounded), so is U. Now the space Y = [— oo, +00] 
is compact, hence has a unique uniform structure, and this is precompact. With 

a multiply connected domain in the complex plane, and $ the set of functions 
M(z, •) (z E 0) above, this gives a uniform structure the completion of which is 
the Martin compactification. The same remark also applies, with Y = P, to the 
family of functions of C >-> F(z,() = d z G(z,()/d z G(a,() for 2 e Q, F(a,a) = 1. 
This gives a uniform structure which we call the gradient structure. Since P is 
compact we have 

Theorem 4.2. The gradient structure Uq is precompact. 

To make the link with the Constantinescu-Cornea theorem we recall that as 
soon as a uniform structure is introduced on a set X, there automatically arises a 
complete uniform space (X,U) and an injection map / : X — > X such that f{X) 
is dense in X and U = f~ x (hi). The triple (/, X,U) is the completion of X with 
respect to 11. Moreover the space X is compact if and only if the uniform structure 
U is precompact. In our setting, with a multiply connected domain f2 c C, we can 
therefore formulate Theorem 14.21 in a more precise way as follows. 

Theorem 4.3. The completion of the space Q with respect to its gradient structure 
Ug yields a compact space OP — il U <9 e f2, the gradient compactification of fl. 

In analogy with the Martin compactification, we call the set d Q £1 = £l e \ ft the 
gradient boundary. In Theorem l3.2l we showed that for multiply connected domains 
with analytic boundary, the gradient boundary is identical with the Euclidian one. 

4.2. An estimate of P. Levy. Thinking of the Martin boundary as a local concept 
we may, in view of Example 13.11 and beginning of subsection 13.31 ask for local 
approximation of the Green's function of a general domain with sufficiently smooth 
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boundary by the Green's function of the upper half-plane. Sharp estimates in this 
direction have been obtained by P. Levy [35]. Below we review Levy's method. 
Let z £ f2, d(z) = dist (z, d£l), c a point on dQ on distance d{z) from z and let z' 

be the symmetric point of z with respect to the boundary d£l, so that c — -(z + z'). 

We assume that z' ^ f2, which is the case for instance if d£l satisfies an exterior 
ball condition (which we henceforth assume) and d(z) is sufficiently small. Let R 
the radius of the largest circle in £1 which is tangent to dil at c and R! the radius 
of the largest circle outside f2 which is tangent to dil at c. We denote by a, a' the 
centers of the respective circles, assuming that R' < oo; the case R' = oo can be 
easily covered by a limit argument. The points a, z, c, z' , a' lie along a straight line. 

Theorem 4.4. In the above notation, 

, ( 2d \ , ,z'-C, ,s , ( 2d 

log i 1 - 2Wd) < log "7^1 - ^ < (i + ^ 

/or every £ G f2. 



Proof. The function £ i— > log | | — G(z, Q is harmonic in Q and equals log | 



on c?S7. The level lines of the latter function are circles, namely Apollonius circles 
with respect to the points z and z' . Recall that the Apollonius circles with centers 
z and z' are r fc = {( G C : \z' - C| = - C|} for fc > 0. 

z' - C 

We conclude from the above that we get upper and lower bounds for log I 1 

z-C 

on dfl by considering one Apollonius circle which lies entirely inside f2 and one 
which lies entirely outside f2. Optimal choices with respect to the given data are 
those Apollonius circles which are tangent to the largest interior and exterior balls. 
The points of tangency are denoted by b and b' . It is easily seen that 



for £ £ Since 



we get, for ( £ dfl, 



|/-6| = 2i? + d, |z-6|=2ii-d, 
|z' - 6'| = 2i?' - d, \z-b'\ = 2R' + d 



lo g 1 ~ 777^^ < log I 7 < log 1 



2R' + dj °'z-C V 2R-d 



and hence 



log ( 1 " 777^) <log|^4|-G(,,C)<logfl ' 



2R' + dJ 01 2- C V 2R-d / 

Now the assertion of the theorem follows from the maximum principle. □ 



5. The Poincare metric, Taylor coefficients and level lines 

5.1. The Poincare metric on level lines for the Green function. Here we 
shall discuss some properties of level lines of Green's function and in particular 
their connections with the Poincare metric. Later on we shall interpret these level 
lines as geodesies for a different metric. 

Consider a simply connected domain f2, let da z — p(z)\dz\ be the Poincare metric 
and G(z, Q the Green's function of f2. 
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Proposition 5.1. For every a G fl, the density p(z) of the Poincare metric is 
(5.1) p{z) - = 



sinh G(z, a) 

In particular, on each level line G(z, a) = c, 

1 dG 
p{z) = -2^-cdn- z M > 

i.e., p is proportional to the harmonic measure with respect to a, or equivalently to 
the Poisson kernel of the enclosed domain. 

Proof. The proof is classical and we recall it for the convenience of the reader. With 
Q(z,a) = G(z,a) +iG*(z,a) an analytic completion of G(z,a) with respect to z, 
the map zh( = e S(z,a) genc j s q on to the unit disk, for which the Poincare metric 
is 

drr ^ 

By conformal invariance 



da z = p(z)\dz\ 



1 _ | e -g( z ,o)|2 



_ e -G(«.a)| 2 flG (z>o) | _ |fg(z,q)| 

I _ e -2G(»,B) - sinh (G(z, a)) 1 U 

which is the desired result. □ 

Remark 5.1. If f2 is multiply connected the above expression for the Poincare metric 
still holds if G(z, £) is interpreted as the Green's function for the universal covering 
surface of f2. This is not single-valued in fl, but the combination appearing in the 
expression for p is single- valued. 

As an application, we reprove a result due to T. Kubo [3B]. With CI simply 
connected as above, let a € O, let p be a positive measure with compact support in Q 
and choose c > so that the support of p is contained in K = {z € CI : G(z, a) > c}. 
We define the Green's potential G M of p by 

G»{z) = J G(z.X)dp(0. 

The functions G, G^ are both harmonic in CI \ K and vanish on dCl. Hence, with 
do~ z the Poincare metric as above, 

/ G»(z)da z = -—?—- [ G»(z)^(z,a)ds 
JdK 2smhcJ dK dn 

= — — / G»(z) — (z,a)ds=—— G (z,a)—ds 

2 sinh c J d (n\K) on 2smhc J din \ K) dn 



[ dG "dc = - p(K) 
JdK 9n sinh c 



2 sinh c <9n 
Thus, under the above assumptions, 

Theorem 5.2. The average of the Green potential G'% with respect to the Poincare 
metric, on a level line G(-,a) — c equals a constant (independent of p) times the 
total mass of p. 
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For the Poincare metric in a simply connected domain we have the estimates 
(see, for example, [4*0] ) 



(5.2) 



1 



< p(z) < 



1 



Ad{z) ~ rx ~' ~ d(z)' 

where the lower bound is a consequence of the Koebe one-quarter theorem. Com- 
pare also l|5.8p below. By (|5.ip this gives the following estimate of the distance to 
the boundary directly in terms of the Green function: 

sinhG(z,a) sinhG(z, a) 
< d(z) < 



I 8G 



(z,a)\ 



-{z,a)\ 



See [IDJ for possible applications of such estimates to computer graphics. 

Also the Bergman kernel can provide estimates for the distance to the boundary, 
even in the multiply connected case. In fact, setting 

nm+n 

K(m ' n) ^0=g^K(z,C) 

we have, according to P. Davis and H. Pollak [10], the Cauchy-Hadamard type 
formula 

(5.3) 



-J- = limsup- (K( n < n \z,z)) 2 " = lim sup ( -^—K {n ^{2 

d(Z) n^oo n \ J n^oo 



for any zeO. 

In the simply connected case, this gives an interesting formula for the distance 
to the boundary. Let <j>(Q = exp(— £/(£, z j) = <zi(C — z) + a,2{( — z) 2 + ■ ■ ■ , a\ > 0, 
be the conformal map from fl to the unit disk which takes a given point z to the 
origin. Then 



(5.4) 



1 



d{z) 



= lim sup 



k=0 



L Ji °»-i+i^jlc=XC) 

j—k 



5.2. Taylor coefficients. Above we saw how the distance to the boundary controls 
the Poincare metric. Below we shall see more generally how this distance controls 
the Taylor coefficients of the Green's function, which in the simply connected case 
embody the Poincare metric. For any multiply connected domain Q, let H (z, be 
the regular part of the Green's function, defined by 

(5.5) G(z,0 = -log\z-(\+H(z,(), 

and let Q(z, C) and 7i{z, Q = H(z : Q) + \H*(z,Q denote analytic completions of 
G(z, C) and H(z,£) with respect to z, for fixed £• Then Q{z,Q is multivalued, 
but Ti(z, C) is a perfectly well-defined analytic function (in z), uniquely determined 
after the normalization Im7i(£, C) = 0, henceforth assumed. Thus Tt(zX) can be 
expanded in a power series around z = (: 



(5.6) 



H(z, C) = c (C) + Cl (C)(z - C) + c 2 (()(z - C) 2 + 



where Co is real by the normalization chosen. 

The first few of the coefficients c n (C) arc domain function which have geometric 
and physical relevance. The constant term, co(C), is sometimes called the Robin 
constant and e~ c °^ is a kind of capacity (if one allows oo 6 D, then co(oo) is the 
ordinary logarithmic capacity of C\ fi), cf. [35]. 

As follows from (|3.ip . cq{z) is related to the Bergman kernel by 

-Aco(z) = 4ttK(z,z). 
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Under conformal mappings cq(z) transforms in such a way that 

da = e- c °M\dz\ 

is a conformally invariant metric (see further Section [7|) . When Q is simply con- 
nected this metric coincides with the Poincare metric and with the Bergman metric, 
normalized to be da = y/irK(z, z)\dz\, but for multiply connected domains none of 
these metrics are the same. Note also that comparison with (|5.1[) shows that 

c (z) = -log ldzK ' l 



' sinh G(z, a) 

for all z, a S f2 when f2 is simply connected. In the multiply connected case we still 
have that the limit, as z — > a, of the right member equals co(a). 

In a simply connected domain any conformally invariant metric has constant 
Gaussian curvature, because the curvature transforms as a scalar and the conformal 
group acts transitively on a simply connected domain. For the above metrics the 
constant curvature is negative: KQauss = — 4. This means that cq satisfies the 
Liouville equation —Acq = 4e~ 2c °, and for the Bergman and Schiffer kernels it 
means that 

OHog K(z,Q „ K( n d 2 log L(z,Q _ J( n 
5.7 — 7 = 2ttK{z, 0, j—r- = -2nL{z, Q. 

We emphasize that these relations only hold for simply connected domains. For 
multiply connected domains there are counterparts of for example (|5.7p involving 
also the zeros of the Bergman kernel, see [5T], [35], [35] and Remark 19.21 in the 
present paper. 

The function cq(z) can be estimated in terms of the distance d(z) = d{z, <9f2) to 
the boundary by 

(5.8) logd(-z) < co(z) < logd(z) + A 

for some constant A. The lower bound is an elementary consequence of the maxi- 
mum principle combined with monotonicity properties of co(z) with respect to the 
domain. The upper bound depends on the nature of the domain. If for example 
the domain is convex, one can take A = log 2, and for a general simply connected 
domain (|5.8[) is the same as (|5.2[) , i.e., A — log 4 works. See [35], pQ for further 
discussions. 

For the higher coefficients we have the following estimates, one of which will be 
used in Section [HJ 

Lemma 5.3. For n > 1, 

Proof. The derivative of the Green's function has the Taylor expansion, with respect 
to z, 

dG(z,() 1 v-^ rs\t rtn-i 

dz = ~7T7 + Z^ ne " (0(^-0 • 

n— 1 

Thus, for any region D C CI containing £, 

1 f <9G(z,C) dz 1 f dG(z,()+i*dG(z,0 



nc (C) = -f 2 dG(Z ' C) dZ — I 



On taking D = {z E fl : G(z, £) > e}, where e > 0, we have dG(-, () = along dD, 
while —*dG(-, () can be considered as a positive measure of total mass 27r on 3D. 
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By letting e — > this gives the desired estimates: 

\nc n {0\ = |— / -7— | < 



^Jan (z-0 n ~ d(C,dQ)- 

□ 

5.3. The Green's function by domain variations. The coefficient c\ in (|5.6[) 

can be directly obtained via the formula 

2 Jan oz 2i ,/ an dn 

which follows from the residue theorem. Instead of the residue theorem one may 
use the fact that 2ci is the gradient of Co, as is easily checked (see (|7.23[) ). combined 
with the Hadamard variational formula, 

If dG(;z) dG(;Q x , 

5Gn{z,Q = — — Snds, 

27r Jan on dn 

which gives the infinitesimal change of the Green function under an infinitesimal de- 
formation Sn of dtt in the outward normal direction. Since SHfi(z, Q = SGn{z, C), 
one obtains <5co(C) by choosing z = £ above, and then the gradient of Co is obtained 
by choosing Sn suitably. See [19] . Lemma 8.4, for further details. 

Here we wish to expand slightly on another use of the Hadamard formula. To 
avoid infinitesimals one may replace 5 by 4r, where t is a time parameter, so that 
% means the velocity of dfl in the normal direction. Of special interest is to take 
this normal velocity proportional to the normal derivative of the Green's function 
itself, say 

Sn _ dG(-,a) 
St dn ' 

where a £ O is a fixed point. This is called Laplacian growth, or Hele-Shaw flow 
with a point source (see, e.g., [55 for further information), and if V(a) = 4z denotes 
the corresponding derivative, acting on domain functionals, the Hadamard formula 
gives 

V(n\r (hA 1 ( dG(-,a) dG(-,b) dG(;c) ^ 

hence that V(a)Go(6, c) is totally symmetric in a, b, c. 

This remarkable fact appears in a series of articles by M. Mineev, P. Wiegmann, 
A. Zabrodin, I. Krichever, A. Marshakov, L. Takhtajan, for example [41], [35], [57] . 
from a more general perspective to be a consequence of two things. The first is that 
the regular part H(z, £) of the Green's function is a double variational derivative 
of an energy functional F = F(£l), which can be identified with the logarithm of a 
certain r-function [34J. Precisely, Ha(a,b) = V(a)V(6)i r (f2), hence 

G n (o, b) = - log \a-b\+ V{a)V(b)F{fL), 

where 

F(Q) = -L f f \og\z- C\dA(z)dA(C), 

87r JD(0,R)\Q JlS(0,R)\(t 

dA denoting area measure and where R is large enough, so that f2 C O(0, R). The 
second ingredient is that the Dirchlet problem is "integrable" in the sense that 

V(a)V(6) = V(6)V(a). 

Clearly these two facts embody the total symmetry of V(a)Gn(b, c). 
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5.4. Level lines of harmonic functions as geodesies and trajectories. Here 
we shall interpret level lines of harmonic functions as geodesies in riemannian mani- 
folds and trajectories of hamiltonian systems. Curvature of level lines and geodesies 
will be discussed in Section [51 

Proposition 5.4. Let u be harmonic in some domain, u* a harmonic conjugate 
of u and let tp > be any smooth function in one real variable. Then, away from 
critical points of u, the level lines of u are geodesies for the metric 

da = ip(u*)\Vu\\dz\. 

Proof. Let $ be a primitive function of tp and let 7 be a level line of u, with 
V11 =/= on 7. By the Cauchy-Riemann equations, this level line is simultaneously 
an integral curve of Vw*, and |Vu| = |Vu*|. Thus, if zq and z\ denote the end 
points of 7 ordered so that u*(zq) < u*(zi), 

da= I <f>'{u*)\Vu*\\dz\ = ( &(u*)du* = - <f>(u*(z )). 

"y J -y J -y 

Since, along a curve in general, |Vu*||dz| > du* , integration of da along any curve 
from zq to z\ will give a value > $(u*(zi)) — <&(w*(z )). Thus 7 is a geodesic. □ 

As a simple remark on trajectories, the level lines of any (smooth) function u 
in a domain fl c C are trajectories of the hamiltonian system with phase space 
f2, symplectic form lu = dx A dy and hamiltonian function u (see [5] for the ter- 
minology). The Hamilton equations then are = = — or, in complex 
form, 

, „ N dz n .du 

^ Tt = - 2l ai- 

Such a kind of hamiltonian system describes for example the motion of a point 
vortex in a plane domain, see [37], [43] and Section [8] below. 

The above equations guarantee that the motion is along the level lines of u. 
Assume now that u is harmonic. Then (|5.9p gives also 

d 2 z d 2 udz d 2 udu dV 

(5 0) dtz~ = ~ d^dt = d^dz' = ~ aF' 

where 

1 .__ ,o n dudu 
V = -2^ = - 2 d-zW 
Notice that V is real-valued and that the right-hand side of (|5.10[) is minus the 
gradient of V. Therefore we can think of (|5.10p as an ordinary newtonian system 
for the motion of a unit point mass in the potential V . Thus 

Proposition 5.5. The level lines of any harmonic function u are, away from crit- 
ical points, trajectories for the newtonian system with potential energy — i|Vw| 2 . 

Remark 5.2. The proposition generalizes to the case that Au = c, c constant, with 
V changed to 

V = -i|Vu| 2 + cu. 



If we want to put Proposition 15.51 into a hamiltonian formulation, the domain 
f2 takes the role of configuration space, while phase space is x C. The kinetic 
energy is T = \ \^\ 2 and the hamiltonian 

H = T + V, 
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to be considered as a function of position q = z and momentum p = 4|. The 
Hamilton equations, written in complex notation, are 

dq dH dp _ dH 

Along a trajectory 1 1— > (p(t),q(t)) in phase space, i? (p, g) is constant, say H = E, 
where E is the total energy. Now, the principle of least action in the form of 
Maurpertuis, Euler, Lagrange and Jacobi (see [3J, [2D]), states that the trace in 
configuration space of a hamiltonian trajectory of constant energy E is a geodesic 
for the Jacobi metric, 

dp = Vfds = y/E - V{z)ds. 

Here ds denotes the ordinary euclidean metric in fi, hence arc- length along tra- 
jectories; in a more general context it would be the metric in configuration space 
induced by the kinetic energy. 

In our case V = — i|Vw| 2 , and since we derived the motion from (|5.9p . T = 
^|Vii| 2 . Thus E — and the Jacobi metric becomes 

dp = ^I Vm I l dz l> 

i.e., the least action principle becomes an instance of Proposition [52] 

The hamiltonian system (|5.11[) or (|5.10[) has of course many more trajectories 
(in configuration space) than the level lines of u. In fact, through any point z there 
is one trajectory starting out with any prescribed speed . (Even if we ignore 
the parametrization, the modulus |4|| of the speed affects the trajectory as a point 
set). Some of these other trajectories (one in each direction) can covered by the 
above analysis by mixing u with its harmonic conjugate, as follows. 
Let u* be a harmonic conjugate of u and let 

w = f(z) = u + iu* , 

so that / is conformal away from critical points of u. Then, for any 9 £ R, |Vw| = 
|2§»| = \f'(z)\ = \e w f'(z)\ = |V(cos(9u - sin0tt*)|, so that V remains unchanged 
if u is replaced by cos 9 u — sin 9 u* . Thus the level lines of all such functions are 
also trajectories, and it is easily seen that they represent all trajectories with total 
energy E — 0. In terms of f(z) we can express the conclusion as follows. 

Proposition 5.6. The trace in configuration space fl of the trajectories on the 
energy surface E = of the hamiltonian system \5.11\) or \5.1U\) are exactly the 
inverse images under the conformal map w = f(z) = u + iu* of the straight lines 
in the w -plane. 

As an application of the above we may take u to be the Green's function of a 
domain Q. C C: u(z) = G(z, a), a € fi fixed. With il' = Sl\{zi(a), . . . , z g (a)}, where 
zi(a), . . . , Zg(a) are the critical points of G(z, a), we conclude from Proposition ^. 41 
that that level sets of G(z, a) are geodesies for 

dai = \VG(z,a)\\dz\, 

and 

|VG(z,a)| ,, 
d(J2= G*(z,a) 

And, by changing the roles of G and G* we see that the level lines of the harmonic 
conjugate G* (z, a) are geodesies for 

G{z,a) 
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(and for da\). In the image region under f(z) — G(z,a) + iG*(z,a), the above 
geodesies for do\ correspond to geodesies of the euclidean metric and the geodesies 
for d<j3 correspond to geodesies for the Poincare metric in the right half-plane. 



6. Doubly connected domains 

6.1. The general domain functions for an annulus. This example is funda- 
mental, it reveals in many respect the essential ideas. Our exposition is based on 
ideas which have been elaborated in [53], and previously in [35]. Since we shall 
only discuss conformally invariant questions it will be enough work with annuli. 
We shall use the notation 

A a . b = {z e C : a < \z\ < b}, 

< a < b, for an annulus centered at the origin. The conformal type is determined 
by the quotient b/a (the modulus), so we can fix either a or b. Note also the 
conformal symmetries of A a ^: z i— > e l6 z for any 9 6 R and z >— > ab/z, the latter 
being the conformal reflection about the symmetry line \z\ — Vab, which exchanges 
the outer and inner boundary components. 

Choosing the annulus to be A% r, where R > 1, we can represent the main domain 
functions in terms of elliptic functions. The multivalued function z t— > t = log z 
lifts the annulus to the strip {t e C : < Ret < logi?}, which then represents the 
universal covering surface of A\,r. The lift map extends to the Schottky double of 
the annulus, which is a torus, and takes it onto the universal covering surface of the 
torus, namely C. The covering transformations on C are generated by t 1— > t+2 log i? 
and t 1 — ► t -\- 27ri, hence we have a period lattice with half-periods 



(6.1) 



logi?, 



One fundamental domain, symmetric around the origin, is 

F = {t e C : - logi? < Ret < logi?, -n < Imt < tt}. 

Here the right half of this corresponds to the annulus Ai,r itself, while the left 
half corresponds to the copy of the annulus which makes up the back-side of the 
Schottky double. The anti-holomorphic involution is the reflection in the imaginary 
axis: J(t) — —i. A homology basis, as chosen in Section [2] can be taken to be 
ot\ = [— logi?, logi?], Pi = [ — irr, i7r] , as point sets. The orientations of a± and (3i 
will however be opposite to the ordinary orientations of the real and imaginary 
axes. 

Sometimes it is advantageous to work in a fundamental domain whose bound- 
ary is made up of preimages of the curves in the homology basis chosen. Such a 
fundamental domain is 

(6.2) Fi = {t S C : < Ret < 21ogi?, < Imt < 2tt}. 

It is often convenient to scale the period lattice so that it is generated by 1 and t, 
with Imr > 0, in place of 2wi, 2uj2- Then r = 0^2/^2, and in our case we have 

We recall the standard elliptic functions associated to the given period lattice. 
The Weierstrass p-function is 

p(z) = p(z; 2wi, 2uj 2 ) = \ + V (7 — - — 2 

Z l ■t-^ \ (Z + LO) 1 UU 1 
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and the £- and u-functions are 

coo = 



1 / 1 1 Z 



°-w=«n( i -^)«p(^+^ 

where in all cases w = 1m\U)\ + 2m,20J2 and summations and products are taken 
over all (mi, ma) e Z x Z\ (0,0). 
The above functions are related by 

p(z) = -C'(z), C(z) = ^l, 

p(z) is doubly periodic, £(z) acquires constants along the periods, 

C(z + 2lu j )=C(z) + 2 Vj , 

and for a(z) we have 

a(z + 2djj) = — a(z)e 2 

(j = 1,2). The constants rjj are given by r)j — ({u>j) and they satisfy the Legendre 
relation 

17T 

771W2 - ?72Wl = — -. 

In our case 771 is real and positive, 772 is purely imaginary and the Legendre relation 
becomes 

m _ m = 1 

logi? in 21ogi?' 
We record also the differential equation satisfied by p(z): 

p'{z) 2 = 4p 3 -g 2 p(z)~g 3 , 

where 

52-60^^, 53 = 140^^. 

As a first step we shall make some of the functions and differentials appearing 
in Section [2] explicit in the annulus case, and then (in a later subsection) we shall 
study the critical points of the Green's function in some detail. When we work 
directly in the annulus, and with the Schottky double of it realized in the same 
plane by reflection, we shall denote points by letters, z, w, a, 6 and similar. When 
we work on the universal covering surface C, with its period lattice, we shall denote 
the corresponding points t, s, u, v. Thus t = logz, s = logw, u = log a, v = logo. 

The two versions of the abelian differentials of the third kind discussed in Sec- 
tion [2] are, on the universal covering surface, given by 

(6.5) v u ^ v (t) = (£(* - u) - C(t -v) + A)dt, 

(6.6) w u _„(t) = (C(t - u) - C(t -v) + B)dt, 

where A and B are constants (depending on u,v) chosen so that w„_„ get 

the desired periods. This means that 

(C(i - u) - ((t - v) + A)dt, ( (C(t - u) - C(t - v) + A)dt 

' — log R J — 7ri 

are both purely imaginary, and 

[C{t -u)- ((t -v) + B)dt = 0. 

log-R 
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By straightforward calculations this gives 

rj\(u — v) Im (u — v) 



(6.7) A 



logi? 2ilogi? 
?7i (u — v) + 7rim 



(6.8) B = 

log it 

where the integer to depends on the location of u, v relative to the preimage of 
a.\ in the period lattice. With u,v £ Fq, then to — if u and v are in the same 
component of Fq \ [— logR, logi?], which is the case, for example, if Im (u — v) = 0. 
Thus we see, in accordance with Lemma 12.11 that A = B when v = J(u). One 
can also achieve to = in more general situations by working in the fundamental 
domain F\, where a.\ is represented by [0, 2 logi?], hence is part of the boundary. 

Next we compute V(t, s; u, v) by integrating v u - v , see (|2.12|) . The calculation is 
straightforward and the result is 

. , ait — u)a(s — v) , 

V(t,s;u,v) = -log U- Lj L 

a{t — V)<J(S — u) 1 

771 Re {it — s){u — v)) Im (t — s)Im (u — v) 
logi? 2 logi? ' 

Pulling this back to the plane of the annulus, i.e., on substituting t = logz etc., we 
get 

v( n i | g(logf)^(logf) | 7?iRe(log^logf) arg ^ arg f 

V(Z ^ aM = - bg l g (logfMlogH) l fog* 2fogl?- 

We record also the final expressions for v a -\, and uj a -b in the plane of the annulus: 



^ , ??ilogf argf dz 



Va-b = (C(log -) - C(l0g T, , ,-, .,. , ,-, • • 

a b logi? 2i log it z 

= (c(log « } ~ c(log b } + fogl? } T- 

The Green's function is by (|2.10p just a special case of V. Choosing s = J(f) = 
— t, w = J(u) = — u and using the first alternative in (|2.10p gives, on the universal 
covering surface, 

. 1, I ait — u)ai— i + u) , rjiKetKeu 

Git, U ) — ~ - lOg ; r ; = . 

v ' ; 2 51 a(t + u)a{-t - u 1 logi? 
Hence, in the annulus (z, a e Ai.r), 

G ( z>0 ) = _ilo g ] ^SfMlogf) , 2^ log|z|log|a| 



! fj(log(za))cr(log^) logi? 

Expanding the u-function in an infinite product gives us what we would have ob- 
tained by the method of images (automorphization) . The result is 

(6.9) G(z,a) = -l g|fc^| - fog | TT g " ~ |>„ |, 

n— 1 v it^* 1 / v za ' 

see g], 0. 

The Bergman kernel is in general 

K(z, QdzdQ = -= dC,. 

it oQ 

Using (|6.6p . (|6.8p . this gives in the present case 

(6.10) K{ t,u) = -(p(t + u) + -\) 

7T logi? 
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for t, u € C, and 
(6.11) 



irza log H 



in the annulus (z, a S Au?). This agrees with expressions derived by Zarankiewicz 
[BD]. See also [IE] (p. 133) and 0]. 

Finally, we elaborate the Schottky-Klein prime function in terms of elliptic func- 
tions (see [5] for representations in terms of Poincare series). It is obtained by 
combining (|2~15| with (|6~6l) and (|6^8)) . The result of that is 

rt 



<j(t — u)a(s — v) 
exp / cj„_„ = — — exp 

a(t — v)a(s — u) 



r]i(t — s)(u — v) 2n\m{t — s) 



\ogR 



\ogR 



which is to be identified with =7 ' — 1 „[ - — \ . Working in the fundamental domain 

E(t,v)E(s,u) 5 

F-y in (|6.2p . which allows m = 0, this gives 



E(t, u) = <j(t — u) exp 



riitu 
tog~R 



/(*)/(«)> 



for some one variable function /. To identify f(t) we differentiate with respect to 
t at t = u: 



dE 
~dt 



\ t=u E(t,u) = </(0)exp 



logi? 



/(«) S 



Since -E(i, u) is to behave like t — uaXt — u the above derivative must be = 1. This 



gives f(u) = ± exp 



Therefore actually 



21ogi? 
i?(i, u) = o{t — u) exp 
and in the plane of the annulus, 

E(z, a) = a(\og -) exp 



2 logi? 



6.2. Eisenstein series. For further need we recall the following classical Eisenstein 
series and other arithmetical functions 1551. 



(6.12) 
(6.13) 



Eo 



2 (r) = 1 - 24]T ^ = 1 - 24 5>(n) g «, 



OO o „ 

n^q 71 



(r) = 1 + 240 73^ = 1 + 240 E CT 3 (n)< 



(6.14) 
where 



E, 



OO ^ n 

rrq n 



(r) = 1 - 504 ]T — ^ = 1 - 504 £ <r 5 (n)<f 



d I n 

and the "nome" q always is related to r by 



q = e 



Thus \q\ < 1. 
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It is well known that Ki(t) and Eq(t) are modular forms of weights 4 and 6 for 
the modular group T = SL2(Z). This is a consequence of Lipschitz' formula, which 
asserts that for any integer k > 2, and with B 2 k being the fc:th Bernoulli number, 

n— 1 

= 1 ~7J^ £ cr fc _i(n)q n = 1 - V ^—^j:- 

B 2k ^ , B 2 k (mT + ny k 

m,n— 1 m.n— 1 ' 

E 2 (t) is not modular but 

nq n 3 
9 



DC- 



(6.15) £5(r) = 1 - 24^ ^- - -^d 3 - 168 9 4 

n=l y ' 

= 1 24q - 72q 2 ~ 96 

7rlmr 

is a non-holomorphic Eisenstein series of weight 2 for T. Actually the series E 2 (t) 
satisfies 

, n .ar + 6, _ , . . ,, 2 6c(ct + d) 

(6.16) E 2 { -) = E 2 {t)(ct + d) 2 + — - - 

CT + d 17T 

for every [ a , J in I\ This is due to the fact that the series 



c 



oo 1 

y t . 

^— ' (mr + n) fe 

n,n— 1 v y 



is not absolutely convergent for k ~ 2. To overcome that difficulty one considers, 
following Hecke, the series, defined for s > 0, 

S| fc (r; S ) = l-^ £ ' 



2fc ("it + n) k \(mT + n\ s ' 

m,n— 1 x 1 



and takes the limit as s — > 0. For fc = 1 this gives (|6.15[) . We will use a similar idea 
in subsection 16.41 to recover the Green's function via an eigenvalue problem. 
We recall the classical discriminant function A, defined as the infinite product 

oo 

A(t) = 9 ]](1- O 24 = q - 24g 2 + 252 9 3 - 1472/ + • • • . 

n=l 

It is a cusp form for T of weight 12 and is related to the Eisenstein series by 

1728A(r) = £ 4 (t) 3 - E 6 (t) 2 . 
Following Ramanujan [48j we also introduce the functions 

oo oo 

(6.17) $rs(q)= ^^m r n s 9 mn . 

Then, 



rn—l n—1 



®rs(q) = U^J $o, S -r(?)< 

We may consider $ r s(q) as a function of r and write $ rs (r). For this we have 
an algebraic relation 



<Mr) =J2 K lmnE 2 (T) l E 4 (T) m E 6 (T) n , 
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for suitable coefficients Ki mn , and where the sum is taken over all positive integers 
I, in, n satisfying 21 + Am + 6n = r + s + 1, I < inf (r, s) + 1. 

It is a fundamental fact |48j that every modular form on T is uniquely expressible 
as a polynomial in £4 and -E^ and that the extension C[£?2, £4, £«] of C[£^, .E^] is 
closed under differentiation, with the following dynamical system of Ramanujan: 

E 2 = 777 ( E % ~ E 4 



(6.18) 



12 
1 



E 4, — ■^(E 2 E i - E 6 ), 
E'q = -(E 2 E 6 - E\), 



where the prime means the differentiation — : — . 

The identity (|6.16p can be expressed by saying that, up to a constant factor, E 2 
defines an affine connection, and the first relation in (|6.18[) similarly says that E4 
is the projective connection which is the curvature of the affine connection E 2 (see 
Section [7| for the terminology) . Hence the Eisenstein series E2 defines a covariant 
derivative sending weight k modular forms / into weight k + 2 modular forms 

f-±v. 

From the system (|6.18[) we obtain 



(6.19) 



E4 — E? — YIE^i 

E 6 =E%- 18E 2 E! 2 + 36E%. 



An important consequence is that the function E 2 is a solution of the Chazy equa- 
tion 

(6.20) E' 2 " = E 2 E' 2 ' - ~E' 2 2 , 

This means that we may use E 2 ,E' 2 ,E' 2 ' instead of E 2 ,E±,E§ to express modular 
forms, that is C[E 2 , E4, E 6 ] — C[E 2 , E' 2 , E 2 \. We shall, in a later subsection, give 
an illustration of this fact to determine the modulus of doubly connected domains. 

6.3. Critical points of the Green's function and zeros of the Bergman 
kernel. The critical points of the Green's function G(z,a) are by (|2.10[) . (|2.11[) 
exactly the zeros of f a _j( a ). For the annulus there is exactly one critical point, 
z = zq(cl), and by symmetry this is located on the same diameter as a. For the 
more detailed investigations one may pass to the universal covering, via t — logz, 
u = log a, -u = J(u) = log J (a). Then, by ([63]) . ([6~7]) . 

(6.21) v u _ J{u) (t) = (((t -u)- C(t + u)+ ^ )dt. 

Thus the equation for the representation t = tc(u) of the critical point is 

(6.22) C(Mu )_ u) _ C(iG(u) + , 2 ) = _^^ j 

log it 

where < Reu < logi?. In the plane of the annulus the corrseponding equation, 
for z — zg(cl), becomes 

fa ~o\ f-n z c(a)^ , n , ,_s 2771 log |a| 

(6.23) C(log ) - C(log z G {a)a) = 



a logi? 

These equations have been analyzed by A. Maria [35] 1 an d the results are summa- 
rized in Theorem 16 . 1 1 below . 
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The zeros of the Bergman kernel can be treated fairly explicitly. By (|6.10[) . the 
zeros of K(t,u) on the universal covering surface are those points t = tx(u) for 
which 

(6.24) p(t K (u)+u) = -t^jz. 

logi? 

Thus 

txiu) — s — u, 

where s solves 

(6-25) p(s) = 

log it 

This equation has two solutions in each period parallelogram. None of them are 
real because p is positive on the real axis, while the right member is negative. In 
fact, it is easy to see that every solution has imaginary part n, modulo multiples of 
2tt. 

Let s = s(R) denote that solution of (|6.25[) which satisfies < Res < logi?, 
Ims = 7r. It depends on R or, equivalently, on r via (|6.3[) . Define also 

(6.26) p = p(R)=e Rcs = \e s \. 

According to [13], s can be explicitly computed to be 

. 3 logi? 3456 V3tt log f io ° $ 5 , 6 (f) - A(t) 

S(R) = ITT + — - — + / — - — —3 (t - T)dt, 

2 5 J T $2,3 W 5 

where the integration is vertically upwards starting at t = lo 'J fl . 

The following theorem summarizes results in |39j and [53j and combines these 
with our findings for the zeros of the Bergman kernel. 

Theorem 6.1. With p defined by H6.26}) we have 

(i) 1< p < VR. 

(ii) The Green's function G(z, a) of A\ # has, for any given a € Ax t u, a unique 
critical point z — zc{a). This is located on the same diameter as a but on 
the opposite side of the hole. More precisely, 

z G {a) = -g{\a\)j- v 
\a\ 

where g : (1, JZ) — > (1,-R) is an increasing function which maps (1,R) onto 
the relatively compact subinterval (p, R/p). It satisfies g{x)g{R/ x) = R for 
1 < x < R, in particular g(y/~R) = \/~R, and 

lim g'(a) — lim g'(a) = 0. 

a— *1 a—*R 

(iii) When p < \a\ < R/p the Bergman kernel K{z,a) has no zeros. 

(iv) When 1 < \a\ < p or R/p < \a\ < R, K{z,a) has exactly one zero, z = 
zk{cl), and this is explicitly given by 



z K (a) 



— ^ when 1 < \a\ < p, 

when R/p < \a\ < R. 



This theorem is special to the annulus but is effective. We do not know of any 
such precise results in higher connectivity. 

Note that the theorem confirms in a precise way the assertion of Corollary 13.31 
that the limiting set for the critical points of the Green's function is exactly the 
zeros of the Bergman kernel when the parameter variable is on the boundary. In 
fact, choosing a on the positive real axis (for simplicity of notation) the theorem 
shows that 

lim za(a) = lim zk(o) = — p, 

a — >1 a — >1 
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R 

lim zg{o) — lim zk{o) = • 

a^B. a^B. p 

In addition we deduce from the theorem the following remarkable dichotomy result. 

Corollary 6.2. The annulus A\n is the disjoint union of the set of critical points 
{zg(o,)} of the Green's function, the set of zeros (z^-(a)} of the Bergman kernel 
and the two circles {\z\ — p} and {\z\ = 

Proof, (of theorem) 

By rotational symmetry we may assume that a is real and positive, namely 
1 < a < R. It has been shown in [39] that 

zg{o) = ~9 (a) 

for some differentiable function g : (1,-R) — > (l,i?), which is increasing, more pre- 
cisely g' > 0, and satisfies 

< lim g(a) < lim g(a) < R, 

a— >1 a— >B 

(the latter follows also from [55] and from our Theorem 13. 1[) and 

lim g' (a) = lim g' (a) = 0. 

a— *1 a— >i? 

Slightly more explicitly we have g(e u ) = e^ u \ where the function / : (0,logi?) — > 
(0,logi?) is defined as the unique solution of 

C(/(u) + m - u) - C(/(«) + iTr + u) = 

or 

(6.27) ic(/(«) + iTr - «) - C(/(«) + i^r + u) = 

lu log It 

Since > 0, and hence /' > 0, /(0) = lim u _>o f(u) exists, and letting u — > in 
(|6.27[) gives (in view of (f = — p) 

^ + ^ = ~^R- 

It follows that 

/(0) + iTT = S (i?) 

(see (|5^5]l - (|5^5[l ). In other words, 

lim g(a) = - e Res ^ = -p(i?). 

a— »1 

Similarly, 

We conclude that p < R/p, i.e., that 1 < p < VR. By this parts (i)and (ii) of 
the theorem are proven. 

The assertions (iii) and (iv), about the Bergman zeros, are easy consequences of 
the analysis made before the statement of the theorem. In fact, since p is an even 
function, the two solutions of (|6.25p are ±s, modulo the period lattice. It follows 
from (|6.10p that, given u £ Fq, the (possible) solutions t — tfc{u) of K(t,u) = 
are represented by 



(6.28) t K {u) = 



s — u, 
— s — u, 



modulo the period lattice. Both u and tx (u) shall correspond to points in A\^r, i.e., 
< Rew < logi? and < Re^(tt) < logi?. This occurs if and only if < Rew < s 
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or log R — s < Re u < log R, and then we have the same inequalities for the zero 
txiu) itself: < Retx(u) < s and logi? — s < Retx (u) < logi?, respectively. 
In the plane of the annulus (|6.28p becomes 



z K {a) 



] cxp(- 



and since this zk{o) is in A\^ if and only if 1 < \a\ < p or Rj p < \a\ < R we have 
now proved statements (iii) and (iv) in the theorem. 

□ 



It follows from (i) and (ii) of Theorem 16.11 that the critical points of the Green's 
function for any annulus of modulus R > 1 are located in the annulus of modulus 
Rj p 2 symmetrically centered in the original annulus. Since p = p(R) > 1 for all 
R > 1 and p depends smoothly on R one easily concludes the following. 

Corollary 6.3. Let R > 1. There exist a sequence {p n }, 1 = Po < Pi < P2 < ■ ■ ■ < 
VR with linin^oo p n = VR such that, on setting A n — {z £ C : p n < \z\ < R/ p n }, 
the critical points for the Green's function of A n are all contained in A n+ \. 

One might ask what would be corresponding statement in higher connectivity. 

6.4. Spectral point of view. The Green's function can also be obtained via spec- 
tral problems for the Laplace operator. This requires a choice of a metric, or at 
least a volume form. For a metric da = p\dz\, the volume form is p 2 dxdy and 

.1 d 2 d 2 

the invariant Laplacian is — A, where A = ——, - + — — . Thus a natural spectral 

p 2 ox 2 ay 2 

problem in a domain is 

{—Am = Xp 2 u in Q., 
u = on dfl. 

If Ai, A2, . . . are the eigenvalues and ui,i«2, . . . the eigenfunctions, normalized 

\u n \ 2 p 2 dxdy = 1, 
then the Green's function of £1 is given formally by 



G(z,w) = 2tt 



E u n (z)u n (w) 
\ 



n=l 

where however the sum may converge only in a weak sense. See for example [5] , [13] . 
(The factor 2ir appears because we have no factor 1/2tt in front of the singularity of 
the Green's function.) One way to cope with the convergence problem is to consider 
the corresponding Dirichlet series 



G s (z,w) = 2irJ2 



u n (z)u n (w) 



which is an analytic function in s for Res > 1. One then studies the behaviour 
as s — > 1. To continue analytically G s (z,w) we shall use the second Kronecker 
limit formula which is basically a connection between Epstein zeta functions and 
Eisenstein series [55] . 

The eigenfunctions u n form an orthonormal set both with respect to the L 2 -inner 
product J n uvp 2 dxdy and with respect to the Dirichlet inner product D(u,v) — 
J n Vw • "S/vdxdy. Note that the A„ and u n depend on p, but G(z,w) does not. 

We shall elaborate the above approach to the Green's function in the annulus 
case, = Ai^, and with metric coming from interpreting A\^ as a cylinder. This 
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metric is equivalent to the natural flat metric on the universal covering surface of 
A\,r. Thus we consider the lift map z i— > t = log z, by which A\ t R gets identified 
with 

C R = {t e C : < Ret < logi?}/27riZ. 

This can be thought of as a cylinder, with boundary 8Cr represented by the vertical 
lines {Ret = 0} and {Ret = logi?} and provided with the euclidean metric da 2 = 
\dt\ 2 . 

On Ai^r itself the metric is given by p(z) = l/\z\ (z £ A\ t r). We shall keep 
some previous notations, for example (|6.3p . and sometimes use r as a parameter in 
place of R. When working in Cr directly, the eigenvalue problem becomes (since 
p = 1 in C R ) 

{—Ait = Xu in Cr, 
u = on 8Cr. 

This problem can be solved by standard separation of variables techniques. The 
eigenvalues then come with two indices, namely 

2 2 

Tfl 7T 2 2 2 2 

(logi?) 2 

for m — 1,2,..., n = 0, 1, 2, Using also negative values of n one can list the 

eigenfunctions as follows (now deviating from a previous notational convention and 
using the variables x, y, z = x + iy on the universal covering surface). 

I miTX 

■sin- -, (m>l, n = 0), 



l^rnn \Zj 



Vtt log R log R ' 

2 mnx 

, 1 D sm l ^cosny, (m > 1, n > 1), 

t/tt log i? logi? 

2 TTITTX 

: sin — sinny, (m > 1, n < — 1). 



> Vtt log i? logi? 

On using exponentials in place of trigonometric functions the above gives, after 
simplifications and still working on the universal covering surface. 

I^mn (z^)Umn (^) 



G s (z,w) = 2n J2 



1 £ 



21ogi? J (n 2 — m 2 T 2 ) s 

1 E 



2 logi? ^ Imr + np 

where z = x + iy, w = u + if. To go on further we shall need Kronecker's second 
limit formula, contained in the following lemma. 

Lemma 6.4. Let r G C, Imr > 0. Let «,u6K\Z be fixed and define for Re s > 1 
the zeta function 



(m,n)eZ 2 \{(0,0)} 



|mr + n| 

T/ien i/ie analytic continuation at s = 1 o/£ T (s;w, u) is 



g27ri(mii+ntj) 



C T ( 1 ; u, v) = 2tt Vim r - 2tt log | ML_HIIl) | 



:S2 
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where the theta function $1 is defined (cf. \2.17\ )) by 



0i(z)=1?i(z;t)=i9 



1/2' 
1/2 



(*;*0 = E 



:3 i7r(n+i) 2 r+27i-i(n+i)(2+i) 



-2]T(-l) n g2(™+3) 2 s i n (2n+ 1) 



7T.: 



= -2gs sinvrz [|(1 - g n )(l - 2 g "cos27rz + q 2n ), 
and where r\ is the Dedekind eta function, 

n>l 

We remark that for ii,»£Z, the function £ r (1; it, u) reduces to the Epstein zeta 
function 

* , s \ - / Imr 

^ S = E | T 12 

(m,n)SZ 2 \{(0,0)} Vl 1 ' 

which is no longer analytic at s = 1 but only meromorphic. Kronecker's first limit 
formula (below) gives the meromorphic continuation (55] . 

Lemma 6.5. The Epstein zeta function Cr(s) has a meromorphic continuation to 
all C with a simple pole at s = 1, given by 

( T (s) = _1_ + 2tt (7 - ln(2\/lmT|77(T)| 2 )) + £)(s - 1) 

where 7 ?s i/ie Euler constant. In particular C T (0) = — 1. 

Continuing the computation of the Green's function we have, using Lemma 16.41 

G(z, w) = lim G s (z, w) 

s—*l 



Imr 
urn 

s-t-l 27T 



E 



(m,n)eZ 2 \{(0,0)} 



g — (a:— «)ImTgin(]/- u) g — (i+u)Im Tgin(j/ — u) 



\mr + n • 



(Imr) 1 8 ( (x-w)Imr y-v. 
1 ^ t( - S; 2^ ' ~2tT' ~ ' 



S->1 27T 



= -log|i?i(- 



|mr + n\ 2s 
[x + u)Imr y — V 



2ir ' 2tt 



(.x — u)Im r (y — t>)r , 



2tt 



2tt 



log |^i ( 



-(.T + u)Imr (y — v)n 



log|M (Z 2 J )T )|+log|^( 



2tt 

(z + u))r, 



2tt 



2?ri " 

= - 10g ^ ( 2l^ )l+10g| ^ ( 21o^ )| - 

This was all on the universal covering. To get back to one simply substitutes 
logz, logw for z, w: 

^/ x , 1 n / log z log iv , 1 , . . ,logz + log it) . 
G(z,w) = -log|i?i(-5- ^)|+l g| 1 ? l (-°- )|. 



21ogi? 

A similar formula is given in [5] . Compare also (|6.9[) 



21ogi? 
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6.5. A remark on the modulus of a doubly connected domain. The connec- 
tion between the preceeding topics and the Ramanujan partial differential system 
Ip8)l or the Chazy equation (|O0|) can be explored in other directions. To give an 
example, we quote a result in [15] showing that the Bergman minimum integrals 
[3] completely determines the modulus of a doubly connected domain. 

In any domain SI, consider the problem of minimizing ||/|| 2 = J n \ f\ 2 dxdy among 
functions satisfying 

//(«) = ••■ = / ( "- 1) («) = o, 

\/W(a)=l, 

where a £ SI is a fixed point and n > 1. Let /„ = f n , a be the unique minimizcr 
and let A„ — X n (a) — ||/ n ,a|| 2 be the minimum value. Then it is easy to see that 
{fn/V~Ki ■ n = 0, 1, ... } forms an orthonormal basis in the Bergman space, and 
hence that 



(6.29) K(z,w) = jr i fn{z \ fn{w) . 

n—0 

Here the dependence on the point a € SI has disappeared, even though each indi- 
vidual /„ and A n depend on a. 

Let da = y/irK^z, z)\dz\ be the Bergman metric of SI. It has Gaussian curvature 
k = KGauss < 0, given in terms of the Bergman kernel and the Bergman minimum 
integrals by 

2 d 2 log K{z,z) \ 2 (z) 
[ ' ttK(z,z) dzdz Ai(z)' 

The last expression follows from (|6.29p (see [15] for details). Set also 

1 d 2 \og{~n{z)) A («)Ai(z) X 2 (z) 



7(z) 



irK(z,z) dzdz ^2(z) Ai(z) 



Both k and 7 are conformally invariant scalar functions. 

Now we specialize to doubly connected domains. Taking in this subsection SI = 
A r ,i (0 < r < 1) as a model, the conformal invariance implies that k(z) and 7(2) 
depend only on \z\ and that they are constant on dA r \ (the same constant on both 
boundary components). These constants, n r — n\gA rl and = j\dA r i-> depend 
on the modulus R. It is shown in 1151 that 



(6.30) 
where 



7r + 2«v + 8 653 - Uag 2 + 120a 3 



(-Kr-4)» ( 52 -12a 2 ) = 



logr 1 

- — . Ot = -—E 2 {T), 
ITT 12 



92 = ^E 4 (r), g 3 = -J^E 6 (r). 



The quantity on the right-hand side of (|6.30[) depends on r or, equivalently, on t. 
Denote it by /(r): 
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It is shown in [T5] that /(t) is monotone as a function of r. Thus /(t) determines the 
modulus. By using (|6.19[) and Chazy equation, we also find the explicit expressions 

, , . ZE 2 + 2E 2 E 2 

Sir) = - 3 » 

TP' 2 

( 6 - 31 ) ; 2 

,// x -9Ef +4E 2 E 2 'E 2 " -5E' 2 3 
J \ T ) = 5 , 

where the prime denotes differentiation — — r — , as before. 

Z7ri dr 

Thus the Eisenstein series E 2 , which is an affine connection, is useful to study 
the modulus of a doubly connected domain. In the next section we will study affine, 
and projective, connections in more generality. 

7. Projective and affine connections 

7.1. Definitions. We shall review some aspects of the theory of affine and projec- 
tive connections. General references here are [51], [22], [23] . 

Let z = f(t), where z, t are complex variables and / holomorphic. We introduce 
the nonlinear differential operators 

, , dz 1 

{M} = log- = -2108-^, 

{ Zj t \ = 4 bg ^ = -2 jf (4==y, 
1 ' s dt 5 dt v j y ^f ! ' 

r , d 2 dz id dz^ 2 f" 3,/", 2 „ r- , 1 ,„ 

The latter is the schwarzian derivative of /. For {z, t}o there is an additive inde- 
terminancy of 27ri, so actually only its real part is completely well-defined. With s 
an intermediate complex variable, so that z = g(s), s = h(t) for some holomorphic 
functions g and h, we have 

(7.1) {z,t} k (dt) k = {z,s} k (ds) k + {s,t} k (dt) k (fc = 0,l,2), 
see [5T]. In addition, one notes that 

{z, t}o = iff z — t + b (/ is a translation), 
{z, t}x = iff z = at + b (/ is an affine map), 

(7.2) {z, t} 2 — iff z= — -t— (/ is a Mobius transformation). 

ct + d 

The above differential operators can also be obtained as limits of derivatives of 
logarithms of difference quotients of the holomorphic function /. To be precise, let 
z, w and t, u be related by z = f(t), w = f(u). Then 

'{z,t} = lim„^ t log^Ef, 
(7-3) l{z,t}i = 2Km a ^ t -§- t logoff, 



i {z,t} 2 = 61im u ^ t ^logf 



-t • 



Alternatively, for the last case, one may introduce a polarized version of the Schwarzian 
derivative by 

2 w — z 

[z,w;t,u] = 6— —log -. 

atau u — t 

Then 

{z,i} 2 = [z,z;t,i\. 



CRITICAL POINTS OF GREEN'S FUNCTION AND GEOMETRIC FUNCTION THEORY 35 



Now, let M be a Riemann surface. An affine structure on M is a choice of 
holomorphic atlas such that all transition functions between coordinates within 
this atlas are affine maps, i.e., such that for any two coordinates z and z for which 
the domains of definition overlap, the relation {z, z}i = holds. This is a quite 
demanding requirement, and not every Riemann surface admits an affine structure. 
Of the compact Riemann surfaces only those of genus one do. 

Less demanding is a projective structure. It is a choice of holomorphic atlas 
on M such that all coordinate changes are Mobius transformations, i.e., such that 
{z, z\i = whenever z and z are coordinates with overlapping domains of defini- 
tion. By the uniformization theorem every Riemann surface admits a projective 
structure. 

Let <p(z)(dz) rn — <fi(z)(dz) m be local expressions for a differential of order m 
on M. Thus the coefficients transform under holomorphic change of coordinates 
according to 4>(z) = 4>(z)(^) m . We shall allow m not only to be an integer, but also 
a half-integer. This requires a consistent choice of signs in the multipliers (jr ) m , 
i.e., requires a choice of a square root of the canonical bundle (which has transition 
functions Such a square root always exists, in fact there are in general several 
inequivalent choices. We refer to [22], [2B], [25] for details about this. Recall 
however (Section 2) that in case M is the Schottky double of a plane domain then 
there is a canonical choice of square root of the canonical bundle, namely obtained 
by choosing the square root of the Schwarz function to be 1/T(z), where T(z) the 
tangent vector of the oriented boundary. This is the choice of square root which 
will be used in the sequel. 

An affine connection on M is an object which is represented by local differentials 
r{z)dz, f(z)dz,. . . (one in each coordinate variable) glued together according to the 
rule 

(7.4) f(z)dz — r(z)dz ~ {z, z}i dz. 

In the presence of an affine connection it is possible to define, for every k S |Z, a 
covariant derivative Vfe from fc:th order differentials to (A; + l):th order differentials 
by 4>(dz) k i-> (Vk4>)(dz) k+1 , where 

(7.5) V fc = ^ - kr<p. 

oz 

The covariance means that if <p(dz) — cf>(dz) k then V k<f>{dz) k+1 = V 'k<f>(dz) k+1 . 

Similarly, a projective connection on M consists of local quadratic differentials 
q(z)(dz) 2 , q(z)(dz) 2 , . . . , glued together according to 

(7.6) q{S)(dz) 2 = q{z){dzf - {z, z} 2 {dzf . 

From (|7.1[) it follows that this law (as well as (|7.4[) ) is associative. In general, we 
do not require r and q to be holomorphic, although our main interest is in the 
holomorphic (or meromorphic) case. 

In addition to the above one may consider also O-connections, quantities defined 
up to multiples of 2-7ri and which transform according to 

(7.7) p{z) =p{z) - {z, z} Q . 

This means exactly that e plyZ ^ is well-defined and transforms as differential of order 
one. 

A projective connection is less powerful than an affine one, but it still allows 
for certain covariant derivatives: for each m = 0, 1, 2, . . . there is, in the presence 
of a projective connection q, a well-defined linear differential operator A m taking 
differentials of order to differentials of order <p(dz)^~ i— > (A m (j))(dz)~ ~ . 
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The first few look, in a local coordinate z, 

AoO) = (f>, 

oz 

. / ,s 9 3 cb n 36 dq , 

A4(0) = ^ + 10g^ + 10^ + (V + 3^)0. 

oz 4 oz:^ az; az oz" 1 

The covariant derivative A m will be called the m:th order Bol operator, for reasons 
to be explained further on. 

Any affine connection r gives rise to a projective connection q by 

dr 1 9 

(7-8) *=T Z -2 r - 

This q is sometimes called the "curvature" of r (cf. [H]). In fact, its definition is 
analogous to that of the curvature form in ordinary differential geometry, see |20j . 
Slightly more generally than (|7.8[) . any two affine connections rj(z)dz, j = 1,2, 
combine into a projective connection by 

l.dn dr 2 . 
q= 2 { ^ + 9z"" rir2) - 
Not every projective connection is the curvature of an affine connection. In the 
case that a projective connection comes from an affine connection, as in (|7.8p . the 
corresponding covariant derivatives are related by 

(7.9) Ai=V , A 2 = ViV_i, A 3 = ViV V_i, etc. 

See [33] for the proof. 

The difference between two projective connections is a quadratic differential, and 
the difference between two affine connections is an ordinary differential. Hence, if 
q(z)(dz) 2 is one projective connection the most general one is q(z)(dz) 2 plus a 
quadratic differential. Similarly for affine connections. 

Now, a central fact is that there is a one-to-one correspondence between holomor- 
phic projective connections and projective structures: given a projective connection, 
represented by a holomorphic function q(z) in a general coordinate z, a projective 
coordinate t is obtained by solving the differential equation 

(7.10) {t,z} 2 = q(z). 

It follows from (|7.1|) , (|7.2p that the set of coordinates obtained in this way are related 
by Mobius transformations. In the other direction, given a projective structure, 
a projective connection is obtained by simply setting q(t) = in any projective 
coordinate t. 

One way to solve (|7.10[) . when q(z) is holomorphic, is to consider the second 
order linear differential equation 

d 2 u 1 , , 

(7.11) _ + _ g(z)u = , 

i.e., K<z{u) = 0, for u considered as a differential of order minus one-half. The 
solutions to (|7.10p are exactly the functions 

t(*) = 

ui(z) 

where iti, u 2 are linearly independent solutions of (|7.11[) . 
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In terms of a projective structure, the meaning of the covariant derivatives A m , 
for holomorphic q, is that in any projective coordinate t corresponding to the con- 
nection, A m simply is the k:th order derivative with respect to t: A m ((j}(t)(dt)^ ::L ) — 
^-q^tP- (dt) ~~ V ~ . The fact that the right member here is covariant under Mobius 
transformation is sometimes called "Bol's lemma" [5J, [33] • The precise statement 
of this lemma is as follows. 

Lemma 7.1. [BJ Let z = f(t) = ad - be = 1, A(t) = (ct + d)- 1 , so that 

f'{t) = A(<) 2 . Then, for any smooth function F(z) and any positive integer m, 

For m = 1 this is the ordinary chain rule, holding for any change of coordinate 
/, whereas for m > 2 the formula holds only for Mobius transformations. 

In any projective coordinate t, a natural fundamental set (basis) of solutions of 
the equation A m u = is {1, t, . . . , i™" 1 }. Considering these functions as ^^j^-forms 

(namely (dt)~ ~ , t(dt)~ ~ etc.) and turning to a general coordinate z, this basis 
transforms into {w™ 1-1 , u" l ~ 2 u 2 , ■ ■ ■ , u™ -1 }, where u\ = *J%, U2 = are ~\~ 

forms. The transformation property of A m then shows that {it™ -1 , u™~ 2 it2, ■ • • , 
is a fundamental solution set for A m u = when considered as a differential equa- 
tion in the z variable. It follows that the operator A m agrees with the (m — l)-fold 
symmetric product of A2: A m = S' m_1 (A2). 

In summary, the Bol operators A m , m > 2, are all generated by A2, and their 
solutions are generated by two solutions of A2U = 0. In the last section of the paper 
we shall discuss a method, via a prepotential, of finding all solutions of A2U = 
from a single one. 

Turning to affine connections there are analogous statements as for projective 
connections: there is a one-to-one correspondence between affine structures and 
holomorphic affine connections, the correspondence between an affine (or "flat") 
coordinate r and an affine connection r(z) given in a general coordinate z being 

(7.12) {r,z} 1 = r(z). 

This equation can be directly integrated as 



t(z) = J exp( / r(()d£)dw. 
One may also think of J7.12[) as 



ViV T 

Ot 



or Vi(dr) = 0, cf. [12]. With r used as coordinate, V m (</)) = for any m- 
differential 0(r)(dr) m . _ 

Along with A m and V m it is natural to introduce their conjugates L m and V m , 
defined by replacing Jj- by ^ and q, r by f. 

Even more powerful than an affine connection is a hermitean metric, which we 
in the present section prefer to write in either of the following two ways: 

(7.13) da=^- = e p M\dz\. 

u{z) 

Here u> > transforms as the coefficient of a form of bidegree (— |,— ^) (i.e., 
uj(dz)~ 1 ^ 2 (dz)~ 1 ^ 2 is invariant), whereas p(z) = — logw(z), in analogy with (|7.4p . 
(|7.6p . transforms as the real part of a 0-connection: 

(7.14) p{z) =p(z)-Re{z,z} . 
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Any hermitean metric gives rise to a, not necessarily holomorphic, affine connec- 
tion by 

(7.15) r = _ 2 |-log W = 2^. 

Oz oz 

This relationship says that the covariant derivative of the metric vanishes: V 1 (1/w) = 
Vi(l/w) = 0, or, which turns out to be the same, 

(7.16) V_iw = V_iw = 0. 
Clearly ([779]) . (f7TT6|) imply that 

(7.17) A 2 uj = 0, 

and also that A20J — 0. Moreover, it is easy to check that A m (jj m_1 = A m ix) m_1 = 
for any m > 2. 

The Gaussian curvature of the metric (|7.13[) is 

(7.18) n Gauss = -e-^Ap = 4(c |^ - |^| 2 ), 

ozoz 1 oz 1 

which is a real and scalar quantity. In view of (|7.15p , it follows that r is holomorphic 
if and only if KQauss = 0. The corresponding projective connection is 

(7-19) »-^ = 2 &-&n=-2*^. 

oz 2 oz Oz uj 

From (HHD, (|7TTg)l follows §f = -2w 2 §§ or, since k is real, 

dn Gauss = ~2uj 2 (^dz + ^-dz). 

oz Oz 

Cf. [23]. It follows that q is holomorphic if and only if the curvature KGauss is 
constant. When this is the case there is a corresponding projective structure, and 
working in any projective coordinate t, equation (|7.17[) and its conjugate can be 
directly integrated to give uj(t) = att + bt + bt + c, a, c real, b complex. The value 
of the curvature comes out to be KQauss = 4(ac — \b\ 2 ). It also follows that if 17 is 
complete for the metric, meaning roughly that uj = on dQ, then in any projective 
coordinate cT2 is a circle or a straight line. 

Being complete and having constant negative curvature characterizes the Poincare 
metric up to a constant factor. Starting from the domain, first the Poincare metric 
(normalized so that KGauss = — 4), then the projective connection q(z) and finally 
projective coordinates, can be obtained directly by solving differential equations 
in CI. In fact, p(z) satisfies by (IT. 18[) a Liouville equation and blows up on the 
boundary: 

(Ap = 4e 2 P in Q, 
1 p = +00 on dfl. 

It can be shown that this boundary value problem, when properly formulated, has a 
unique solution p (see, e.g., PQ, Ch.l). From this solution, q is obtained via (|7.19[) . 
and finally projective coordinates t are gotten by solving (|7.10|) . more precisely 
by taking the quotient between two solutions of (|7.11|) . The above indicates one 
of the early attempts to solve the uniformization theorem. It was proposed by 
H. A. Schwarz and later brought to an end by Picard, Poincare and Bieberbach; 
see [5]. 

This approach has recently gained new interest in some versions of string theory, 
in which the above Schwarzian derivative q(z) — {t, z} 2 of the lift map t to the 
universal covering surface takes the role of being the energy-momentum tensor, 
which hence is a projective connection. See for example [IS], [55]. 
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7.2. Examples of projective connections. If £1 is a plane domain bounded by 
finitely many analytic curves, several projective structures and connections can be 
naturally associated to f2. 

Example 7.1. The trivial one, q(z) — 0, i.e., with the coordinate variable z in C as 
a projective coordinate. 

Example 7.2. Let t : Q — > f2 c i rc be a conformal map of ft onto a domain f2 c i rc 
bounded by circles. This map is not uniquely determined, but any two such maps 
are related by a Mobius transformation. Hence it defines a unique projective struc- 
ture for which t is a projective coordinate. This projective structure on Q, extends 
to a projective structure on the Schottky double Ct, as will be seen more exactly in 
the next section. The associated connection coefficient is obtained from (|7.10p . 

Example 7.3. The universal covering surface of Q is conformally equivalent to the 
unit disk. Let tt : D — > Q be a universal covering map. The inverse of tt is 
multivalued, unless fl is simply connected, but the different branches of tt" 1 are 
related by Mobius transformations. Hence a unique projective structure on Q is 
obtained by using local branches of tt" 1 as projective coordinates. The connection 
coefficients q(z) are related to the (branches of the) multivalued liftings tt" 1 by 
(jTrUj) , that is 

{k~ 1 {z),z} 2 = q(z). 

Example 7.4. Let Cl be the Schottky double of f2 and let tt : U — > fl be a universal 
covering map for O. This U can be taken to be either the Riemann sphere (if 
f2 is simply connected), the complex plane (if fl is doubly connected) or the unit 
disk (if f2 has at least three boundary components). In any case, Q, and hence 
fl, is provided with a unique projective structure by using local inverses of the 
uniformization map as projective coordinates. If Q is doubly connected, it even 
gets an affine structure. 

In general, the projective structures in the above examples are all different. For 
example, if Q is the annulus ^4i,_r, then by straight-forward computations one finds 
that the projective connections are given by 

a 

? z = -o, 

z l 

where a = in the cases 1) and 2) above, a = j in case 3) and a = j(l+ n^mi ) m 
case 4). If VI is replaced by a noncircular doubly connected domain also the cases 
1) and 2) will be unequal. 

Example 7.5. Connections also come up from Taylor coefficients of regular parts 
of certain domain functions. For example, for the Taylor coefficients of the regular 
part of the complex Green's function, see Section HT21 we have the following. 

Proposition 7.2. In the notation of Section \5. 6 A plus \3.J$ , the quantities 

(7.20) p(Q = -00(C), 

(7.21) r(C) = -2 C1 (C) = 2^, 

(7.22) q(C) = ~6(^P - 2c 2 (C)) - fm£(C, () 

transform under conformal mapping as, respectively, the real part of a ^-connection, 
an affine connection and a projective connection. 
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Proof. The Green's function transforms under conformal mappings / : f2 — ► Q, 
z = ,f(z), as G(z, Q = G(z, Q, hence 

n(z,o = - io g ^ + n(z,0- 

Now the assertions follow by easy computations, using the formulas (|7.3[) : 

(7 9* (t\ o g g(z,Q | 9 Tiff ^ 8 (f\ 

(7.23) Cl (o = —q^U = 2 ^H*=c = ac H «>0 = dt c °^ 

aCl(0 2c 2 ( C ) = ^U = 2™l =-^( C ,C). 



d( dzd( lz =< dzdC, lz =< 



□ 



We conclude from (17^21 that da = e- C0 ^\d(\ is a conformally invariant metric 
and that r(z) is the associated affine connection. However, q(z) in (|7.22[) is in 
general not the same as the projective connection associated to r(z) by the general 
receipt (|7.8[) . namely 

Q(0^4kC) 2 = -2^P- 2ci (C) 2 . 

In the multiply connected case Q(Q is not holomorphic, whereas q(() in (|7.22p is 
always holomorphic. 

Example 7.6. Coefficients of linear differential equation in f2 transform in compli- 
cated ways under changes of coordinates, and in some cases exactly as connections. 
For example, it is well-known that a second order equation always can be written 
on the form (|7.1ip . that is u" + ^Q(z)u = 0. Then Q(z) works as a projective 
connection, and the corresponding projective structure uniformizes the equation. 
This differential equation will be further discussed in Section Q21 

Example 7.7. Let us spell out the relevant quantities for the unit disk and upper 
half-plane provided with the Poincare metric: 

(i) Unit disk: 

2z 

uj(z) = 1 - \z\ 2 , r(z) = j-Tj, q{z) = 0, Kcauss = -4. 

1 - \z\ 2 

(ii) Upper half-plane: 

oj(z) = 2y, r(z) = -, q{z) = 0, K G auss = -4. 

y 

Note that r(z) is singular on the boundary, while q{z) is not. 

8. Connections on the Schottky double of a plane domain 

8.1. Generalities. If the Riemann surface is the Schottky double of a plane do- 
main f2, then the connection coefficients can be described in terms of pairs of 
functions on ft, in analogy with the previous description (|2.5I) of differentials of 
any half-integer order. Recall (|2.1[) that the transition function between the coor- 
dinate z = z on the back-side f2 and z on the front side f2 is given by the Schwarz 
function. Writing the unit tangent vector as T(z) = jr = e lfl \ the curvature of dfl 
is k = ^| = —iT'(z), where the prime denotes the complex derivative J-j for the 
analytically extended T(z). Then, along dfl, 

{S{z),z}o = -2logT(z) = -2iO, 
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T'(z) 

{S(z),z}idz = -2— ±-!-dz = -2\d9 = -2inds, 
1 (z) 

{S(z), z} 2 dz 2 = -2^^-dz 2 = -2i^cfa 2 . 

T{z) ds 

Therefore, a O-connection on Q can be described as a pair of functions, Pi(z) (for 
the front side) and p 2 (z) (for the backside), on CI, satisfying on dCl the matching 
condition 



Pi(z)=p 2 (z)-2i6. 

Similarly, an affine connection is represented by a pair r 2 (z) satisfying the 

matching condition T\(z)dz — r 2 (z)dz — 2inds, or 



(8.1) r l {z)T{z) = r 2 {z)T{z)-2in onSfl, 

and a projective connection by a pair qi(z), q 2 (z) with 

. dn 



qi(z)T(zY = q 2 (z)T(z) 2 - 2i— on dCl. 

ds 

Note that if z and z are both projective coordinates then q\ = q 2 = and it follows 
that 4^ = on dCl, i.e., that is a circular domain. 

We may notice that (|8.1[) is consistent with the general fact [5T| that the sum 
of residues of a meromorphic affine connection on a compact Riemann surface of 
genus g equals 2(g — 1) (our definition of connection differs from that in [51j by a 
minus sign). In fact, in the case of the Schottky double the sum is, by (|8.1[) . 

-^r / rx{z)dz + -^—. [ r 2 (z)dz = -— [ nds = 2(g-l) 
27Tl Jan 27T1 J -an * Jon 

since there is one outer component and g inner ones. 

Naturally, one may be particularly interested in symmetric, or "real" , connec- 
tions, namely those for which pi = p 2l r\ = r 2 or q\ = q 2 (respectively). The above 
matching conditions then become, for the single representatives p(z), r(z), q(z) on 
CI, 

lmp(z) = —6 on dfl, 
(8.2) Im (r(z)T(z)) = -k on dtt, 

lm(q{z)T(zf) = on dfl. 

ds 

8.2. Formulas for the curvature of a curve. The differential parameter {z, w}i 
turns out to be a useful tool for summarizing various formulas for the curvature of 
a curve in the complex plane. 

The definition of {z,w}\, with z = f(w) analytic, can be written in differential 
form as 

dz 

{z, w}idw — dlog —— — dtlog dz — log dw). 
dw 

Hence 

Im ({z, w}idw) = d&vgdz — dargdw — *rf(log \dz\ — log Idwl). 

To interpret this formula one should let w run along a curve T, say with parametriza- 
tion t i — ► w(t), f£l. Then z(t) = f(w(t)) runs along f(T) and one just replaces 
dz and dw by z'(t) and w'(t), respectively. 

Denoting the curvature of T by k w and the curvature of f(T) by k z we thus have 
the following formula. 
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Proposition 8.1. Under a conformal mapping f : w, the curvature of a curve 
r and its image curve f(T) are related by 

(8.3) Im ({z, w}idw) = n z ds z — n w ds w . 

Here ds z — \dz\, ds w = \dw\. 

Example 8.1. (The curvature in terms of a real parameter.) With T — R and 
z = f(w), so that z = f(t) with t = Rew parametrizes f(T), we have n w = 0, 
which gives the well-known formula 

, d dz dw . ds w 
k z = Im( — (log— ) — ) + k w — 

dw dw ds z ds z 

f(w) dw . 1 f{t) 

m V'H \f'W\\dw\ ) \f(t)\ m /'(*)• 

Example 8.2. (The curvature in terms of an angular parameter; essentially Study's 
formula [46] , p. 125.) With T = oB, w = e 1 *, dw = iwdt we have, 

dsu, ,f{w) iwdt 1 

k 2 = Im({z,w}i — ) + — =Im(- 



ds z > w ds z y f'(w) \f'(w)\\dt\' \f'(w)\ 

1 . wf"(w) 

-(Re J \ J +1). 



I/'HP /'H 

Example 8.3. (The curvature of a curve given as a level line of a harmonic function.) 
Let r be any level line of a harmonic function u in the z-plane. Assume that the 
curve is nonsingular, so that Vu =/= on T. With v any harmonic conjugate of u, 
w = u + iv = g(z) maps T into a vertical line in the w-plane. Thus k w = and the 
curvature of T is obtained from 

K z ds z = Im ({z, w}idw) + n w ds w = — Im ({w, z\\dz) 

, d 9w, , . ,„ d . ,, n du.. , . 

- Im ^ (l0g ^^ ) - Im(2 ^ 1 ° g(|2 ^ l) ^ 

9 9 9 

= -tt log |Vu|dy + — log |Vu|da; = -— log |Vu|ds z . 
ax ay an 

We may may also extract from the above the formula (cf. 3 lj ) 

9V-Y 

In summary, for a level curve T of a harmonic function u the curvature is given 

by 

,„ d . ,_ , d , .du. d 2 u ,du 

8.4 ^ = -— log V« - -— log I — | = —r-i -T' 

an dn dn dn z dn 

In the last expressions we used in place of ^ because the latter becomes obscure 
when applied twice (^ is not really a partial derivative since n is not a coordinate 
of a coordinate system; -4- should be interpreted as a directional derivative along 
the straight line in the normal direction). 

In higher dimensions, (|8.4j) gives the mean curvature of a hypersurface given as 
the level surface of a harmonic function. This can easily be deduced from the more 
general (and well-known, cf. [20) ) formula K mcan = div -S^j, for the mean curvature 
of a level surface of any smooth function u. 
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8.3. Geodesies. The equation for geodesic curves z — z(t) for the metric (|7.13[) is 

where r(z) is the corresponding afhne connection (|7.15p and the parameter t mea- 
sures arc-length with respect to the metric. This equation is the usual geodesic 
equation in differential geometry |20j . just written in complex analytic language. 
The classical Christoffel symbols Tj k (i,j,k = 1,2) turn out to coincide with the 
components ±Rer(z),±Imr(z). An intuitive direct derivation of (|8.5p goes as fol- 
lows. 

The tangent vector along the curve is 4| and the geodesic equation is supposed 
to say that this propagates by parallel transport, i.e., has covariant derivative zero 
along the curve. Since a vector in one complex variable can be thought of as 
a differential of order minus one the relevant covariant derivative (J73J) will be 
V_i = + r(z). This is the covariant version of J^, and the covariant version of 
j- t then is Jj + r(z)^f . Applying this to ^§ gives (|53|) . 

It is convenient to write (|8.5p on the form 

d , dz , s dz 

(8 - 6) dt l0 ^ + riz) Tt-° 

and to decompose it into real and imaginary parts. The real part just contains 
internal information about how the curve is parametrized, namely saying that t 
measures arc-length with respect to the metric. The information about the ge- 
ometry of the curve (that it is a geodesic) is entirely contained in the imaginary 
part, 

(8.7) | arg | + Im(r(2) | )= 0. 

In equation (|8.7p . t can be taken to be any parameter, for example euclidean 
arc-length s. Then the first term equals the euclidean curvature of the geodesic 
and is the unit tangent vector. Thus we have 

Proposition 8.2. The curvature k — k(z, of the geodesic passing through a 
point z and having direction (a unit vector) is given by 

(8.8) «(,, £) = -MK*)£) = -Im (2§^) = * 

ds ds dz ds on 

where ^ denotes the derivative in the rightward normal direction of the curve. In 
particular, the sharp bound 

(8.9) W*,%)\<\r{z)\ 

holds for all geodesies passing through z, and equality is attained if and only ^ is 
tangent to the level line of p at z. 

The above estimate (I8.9[) can be combined with geometric estimates of r(z). For 
example, for the coefficients of the Green's function, as in Section [5] and Exam- 
ple [731 we have by (|7.20p . (|7.2ip that r(z) — —2c\{z) is the affine connection for 
the metric da — e~ c °\dz\. For c\{z) we have the estimate in Lemma 15.31 so the 
above proposition shows that 

|r(z)l " d&dQ) 

for this connection. When is simply connected the metric in question is the 
Poincare metric. Combining with (|8.9p we therefore have the following. 
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Corollary 8.3. For a simply connected domain provided with its Poincare met- 
ric, the curvature for any geodesies through a point z is subject to the estimate 

dz 2 

lK[z 'Ts )l - d(z,dny 

It is allowed here that Q C P contains the point of infinity, and at least in this 
generality the corollary is sharp. Indeed, with £1 = P\B(0, e), e > small, the circle 
with center at any finite point c € f2, and having radius |c|, is almost a geodesic for 
the Poincare metric in f2. The curvature for this circle is l/|c|, hence the bound in 
the corollary is essentially attained at the point 2c on the circle. 



8.4. Some applications in physics. Projective and affine connections come up 
naturally in both classical and modern physics. As to modern physics, e.g., con- 
formal field theory and string theory, we have already mentioned the example with 
the energy-momentum tensor as a projective connection. For further examples, see 
|12j . [58] and references therein. 

For classical physics, we shall briefly mention one example from vortex dynamics. 
Consider in fl an incompressible inviscid fluid which is irrotational except for a point 
vortex of unit strength at a point zq. If Q is simply connected this makes the flow 
uniquely determined. The stream function ip will, at each instant of time, coincide 
with the Green's function, ip(z) — G(z, Zq). However, the flow will not be stationary 
because the vortex will move with the speed obtained by subtracting off, from the 
general flow, the rotationally symmetric singular part corresponding to — log \z— Zq\ 
in the stream function. Thus in fact ip — tp(z,t) = G(z, Zo(t)) and one deduces, in 
the notation of Example I7.5I that the vortex moves along a level line of co(z) and 
that the velocity is given by 

(8-10) ^ = -i^)) = -i^(,oW). 

This means that cq(z) is a kind of stream function for the vortex motion, called 
the Routh stream function, [37], [43] . It also follows that the vortex motion is a 
hamiltonian motion, determined by the symplectic form dx A dy (area 2-form) and 
hamiltonian function H = ^Cq. Note that (|8 . 1 0[) is a special case of (|5.9I) . 

If Q is multiply connected one need to prescribe the circulations jj around the 
holes to make the flow uniquely determined. These circulations will be preserved in 
time (Kelvin's theorem) and everything will be as in the simply connected case but 
with the ordinary Green's function replaced by the hydrodynamic Green's function 
G 7 (z,zo), to be discussed in Section[5] 

A related physical application comes from electrostatics. Think of the comple- 
ment K = C \ f2 as a perfect conductor and let it be grounded, so that its potential 
is zero. Then consider a unit charge located at zq € J7. This will induce charges 
of the opposite sign on dK, namely distributed so that the density with respect 
to arc- length is given by the normal derivative dG g^ Zo ^ ■ These charges will exert a 
force on the charge at z$, and this force is 



F{z ) = -—c 1 {z {i)). 

If Q, is multiply connected and the components of K are not grounded to a common 
zero, but are isolated from each other, then the hydrodynamic Green function shall 
be used in place of the ordinary one, with the jj proportional to the respective 
total charges isolated on the different components of K. 
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9. On Neumann functions and the hydrodynamic Green's function 

9.1. Definitions. Recall that the Bergman kernel and the Schiffer kernel are given 
in terms of the ordinary Green's function by (|3.1|) , (13. 2|) . Decomposing the Green's 
function as in (|5.5p we also have, for the "£- kernel" (|3.4[) . 

7r azoQ 

The reduced Bergman kernel (see Proposition 19. 2\i and its adjoint have several 
similar representations in terms of Neumann type functions and hydrodynamic 
Green's functions. Here we shall briefly review these matters. See [33], [HI] for 
more details. 

By a Neumann function we mean a domain function N a (z, () with a logarithmic 
singularity at a given point £ € fl and satisfying Neumann boundary data given by 
a boundary function a subject to 



(9.1) / ads = 2tt. 

Jan 

The requirements on N a (z, Q are, more exactly, 

(9.2) N a (z, C) = - log l-z - CI + harmonic in fl, 

(9.3) = a on oil, 

on 



(9.4) / N a (-,C)ads = Q. 

Jan 



The final condition (19.41 ) can also be written 

N a (;()*dN a (;z) = 0, 

/an 

and is a normalization which guarantees that iV Q is symmetric: 

N a (z,() = N a (t,z). 

A hydrodynamic Green's function (or "modified Green's function") [37], [19] . 
[9] is defined in terms of g + 1 prescribed circulations, 70, ■ ■ ■ 7 g subject to the 
consistency condition 

(9.5) 7o + ---+7 g = 27r. 

Let 7 = (70, . . -7 g ) denote the entire vector of periods. The defining properties of 
the hydrodynamic Green function, denoted G7 7 (z,£), are 

(9.6) G-y(z, C) = — log \z — C| + harmonic in fi, 

(9.7) G 7 (*,0=&i(O forzer,, 

(9.8) - jf ds(z) = 7j (j = 0, . . . , g), 

(9.9) E7,7MO=0. 

3=0 

Here &_,-(£) denote "floating constants" (they cannot be preassigned), i.e., (|9.7[) 
really means 

dG 7 (-,C) = along 90. 
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Condition (19.91) is a normalization which can be written 

G 1 (;0*dG 1 (-,z) = Q 

I an 

and which guarantees the symmetry, 

G 7 (*,0 = G 7 (C,z). 

The hydrodynamic Green's function can be constructed from the ordinary Green's 
function by 

g 

G 7 0,C) = G(z,C) + ^2 c kj u k(z)uj(Q, 

k,j=0 

where Uk, k = 0,...,g are the harmonic measures. With the above 'Ansatz', 
the requirements (|9.6p and (|9.7p are automatically satisfied, and (|9.8p , (|9.9p give 
a system of equations which determine the coefficients Ckj uniquely. The matrix 
(ckj) is symmetric and positive semidefinite. 

9.2. Reproducing kernels for Dirichlet and Bergman spaces. The Neumann 
and (hydrodynamic) Green's functions have logarithmic singularities, hence cannot 
themselves be reproducing kernels for any Hilbert spaces of harmonic functions. 
However, the singularities disappear when subtracting them, and also after appli- 
cation of the the differential operator . In these cases we do obtain reproducing 
kernels for important spaces. Below we elaborate on these matters, slightly extend- 
ing the analysis in [1] and |51j . 

Let D(u, v) denote the Dirichlet inner product: 



du A* dv, 



D(u,v) = I Vu-\Jvdxdy = I 
Jn Jn 

let a be boundary data as above, satisfying (|9.ip . and define the period vector 7 by 
(9.10) 7fc = / ads. 



Then (|9.5p holds, by (|9.ip . Define H(£l) to be the Hilbert space of all harmonic 
functions u in f2 which satisfy D(u,u) < 00 and 



(9.11) / uads = 0, 

and let H e (Q) be the subspace consisting of those functions which in addition satisfy 

(9.12) l^ dS = ° (J = °" --' S) - 

Note that (|9.1ip just fixes the additive constant in u which the inner product leaves 
free. Except for this additive constant H(fl) does not depend on a. The following 
is a slight extension of results in Ch.V:3 of [3]. 

Proposition 9.1. The reproducing kernels for H(£l) and H e (fl) are, respectively, 

Z7T 

1 

2~TT 

In other words, fc(-,C) € H(Cl) and 

(9.13) «(C) = D(«,fc(-,0) (ueH(n)) 
and similarly for k e (-,£). 



k e (z,0 = —(N a (z,0-G^z,Q). 
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Proof. All verifications are straightforward. Let us just show, for example, that 
(|9.13p holds for fc e (-,C) an d u £ H e (£l)). We may assume that u is smooth up to 
the boundary. The proof amounts to standard applications of Green's formula to 
functions with a singularity. For the Neumann function we get 

D(u,N a (-,0) = f duA*dN a (-X) 
Jn 

= lim / d(u A* dN a (-, 0) = — / uads + 2wu(() = 2nu(C), 



Q\D(C,e) JdQ. 



and for the hydrodynamic Green's function, 

D(u, G 7 (-, 0) = D(G 7 (;{), u) = lim / d(G 7 (-, C) A *du) 



f2\D(C,e) 



= I>(0 / ^ds-Km [ G 7 (- 1 ()pds = 0. 

Now dnH3> follows. □ 

By definition of H e (il), all functions u £ H e (Q) have single- valued harmonic 
conjugates. Thus we can form the space of analytic functions / = u + iu* with 
u £ H e (Q) and with also u* normalized by (|9.11|) . This is simply the ordinary 
Dirichlet space with normalization (|9.11[) . i.e., the Hilbert space A(il) of analytic 
functions in Q provided with the hermitean inner product 

{f,g)-i= [ f'g'dxdy = -— [ dfAdg, 
Jn Zl Jn 

and subject to the normalization 

(9.14) / fads = 0. 

Jan 

The reason for indexing the inner product by —1 is that A.(Cl) shortly will be iden- 
tified as one space in a sequence of weighted Bergman spaces with inner products 
in general denoted (f,g) m , m£*L. We note that, with u = Re/, v = Reg, 

D(u,v) = f du A * dv = -Re ^ f df A dg = Re (/, g)-\. 
Jn 2i Jn 

The reproducing kernel for _A(f2) is the analytic completion fC{z,Q of k e (z, 0, 

normalized by (|9.14p . In terms of the multivalued analytic completions jV a (z,0i 

Gj(z, C) °f N a (z, C) and G 1 (z, £) we therefore have 

1 

27r' 

An adjoint kernel may be introduced as 



^,0 = ^(^,0-^0)- 



c(zX) = ^(M a (z, o+a 7 («,0). 

Z7T 

Remark 9.1. A/" a (-Z,C), £7(2, are n °t symmetric with respect to z, C and are not 
analytic or antianalytic with respect to Q, However, K,(z, is (of course) hermitean 
symmetric and is antianalytic in £. The adjoint kernel, C(z, £) is multivalued ana- 
lytic in both z and £■ See [5T] for more details. 

Example 9.1. For the unit disk, £1 = ED, with a = 1 and (necessarily) 7 = (70) = (2ir) 
we have 

N a (z,() = -\og\z-(\-log\l-zt\, 
G 7 (2,C) = -log|z-Cl+log|l-zC|, 
k(z,() = --log\l-z{\, 
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/C(z,C) = --log(l-zC), 

7T 

£(z,Q = -~log(z-Q 

TT 

N a (z, = ~ l0g(z - C) - l0g(l " ZO, 

e 7 (2,o = -iog(^-0 + iog(i-^0, 

Of course, G 7 (z, () = G(z, (), G-y(z, () = Q(z, () in the simply connected case. 

We shall denote the ordinary Bergman space by B(fl). It is the the Hilbert space 
of analytic functions in Q with (/, /)i < oo, where 

(f,d)i = / fgdxdy 
Jn 

is the inner product. The exact Bergman space, B e (Q), is the subspace consisting 
of those functions of the form F', where F is (single- valued) analytic in fl. The 
Dirichlet space is related to the exact Bergman space by differentiation: A.% : f t—> f 
(or / i— > df) is an isometric isomorphism Ai : A(fl) — > 5 e (fi). 

Proposition 9.2. The reduced Bergman kernel, i.e., the reproducing kernel for 
B e (Q), is given by 

= 2 d 2 N a (z,Q = 2 d 2 G 1 (z,() = d 2 JC(z,Q 
' tt dzd( tt dzdC dzdC 

for any choices of a and 7 as above. Similarly, the corresponding adjoint kernel is 

2 d 2 N a (z,() = 2 d 2 G 1 (z,Q = d 2 C(z,Q 
elz ' U tt dzd( tt dzd( dzd( ' 

Proof. The proof is essentially well-known and straight-forward. Cf. [4] , [52] ■ □ 

Remark 9.2. A beautiful argument, due to M. Schiffer [50 , shows that the function 

F(z) — J L e (t, Qdt is univalent and maps ft onto a domain DcP with 00 S D 

and such that each component of P \ D is convex. See |491 for further information 
and several related issues. 

A particular consequence is that L e (z,Q has no zeros in f2, hence that K e (z,Q 
has 2g zeros, Wi(C)> ■ ■ ■ , ^2g(C); m ^- I n terms of these the following generalization 
|51j . [55] of (|5.7p to the multiply connected case holds 

dzdC 

s fe=i 
For the full Bergman kernel K(z,£) there are similar results [59], but they are 
slightly more complicated because the distribution of the 2g zeros between K(z, () 
and Liz, C) is less clear in this case. For example, Theorem 16.11 shows that the 
number of zeros of K(z, £) in general depends on the location of ( 6 Q. 

9.3. Behavior of Neumann function under conformal mapping. To describe 
how an object (e.g., the Green's function) on a domain Q transforms under confor- 
mal mapping is equivalent to telling what kind of object it is (function, differential, 
connection etc.) when the domain is considered as a Riemann surface. In fact, 
a conformal map onto another domain can be considered simply as a change of 
holomorphic coordinate on the Riemann surface. Many objects associated to a 
domain actually extend to the Schottky double, and since this Riemann surface is 
described by an atlas with two charts overlapping on the boundary of the domain, 
complete information of the conformal behavior of the object then is contained in 
the transition formula on the boundary. 



(9.15) ^— log K e (z, C) = 2wK e (z, Q+nJ^Le (z, m fc (C)) m' fc (C). 
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The ordinary Green's function G(z, £) extends directly to the Schottky double 
as an odd function in each variable, see (|2.10[) . and its differential is the real part 
of normalized abelian differentials of the third kind on the double. Precisely, by 
([230]) . (j2TTj) . Lemma |2~T1 



(9.16) dG(-,Q = -Re<u c _j (f ) = -Rew c _, /(c) . 

When trying to extend the hydrodynamic Green function as an odd function on 
the double, there appear jumps (=2bj) on 9f2. However, the differential <iG 7 (-,C) 
extends perfectly well, along with the conjugate differential *dG 1 (-,()- Therefore 
dG 7 (-,() + i*dG 1 (■,£) is an abelian differential of the third kind with poles at C 
and J(C). Computing the periods gives 

(dG 7 (z,C) +i*dG 7 (*,0) = 2(6,(0 - 6o(0). 
(dG 7 (z,Q + i*dG- f (z,()) = i7j, 

ft- 

i = i, ■ ■ • ,g- 

If we in particular choose the period vector to be 

7= (2vr,0,...,0) 

we see that the /^-periods (j = 1, . . . , g) of dG 1 {-, C) +i *dG 7 (-, C) vanish. Therefore 
(in this case) 

dG 7 (-,C) + i*dG 7 0,C) = -^_ JCC) , 

where dv-b in general denotes the abelian differential of the third kind with the 
same singularities as w a -b (and v a ^b), but normalized so that the /3j-periods vanish. 
Thus, in analogy with (|9.16p we have 

dG 7 (-,C) = - Re ^C-J(C) 

for 7 as above. 

For the Neumann function the situation is slightly more complicated than for the 
Green's functions. We shall consider N a {z, C) and N a {z, C) = N a (z, () + \N*{z, C) 
together with the meromorphic differential 

r a (z, Qdz = 2 dz = dz. 

az oz 

Along the boundary (with respect to z), 

lxa.{T a {z,Qdz) =lm(dAf a (z,Qdz) =* dN a (z,() = -ads. 

In order to discuss questions of conformal invariance the boundary function a 
has to be linked to the domain f2. The most naive choice, a — constant, turns 
out not to give any good behavior under conformal mapping. More promising is 
to relate a to the curvature k of the boundary. If we simply take a = k then (|8.2[) 
shows that T a (z, Cjdz extends to the Schottky double as a symmetric meromorphic 
afBnc connection. Moreover, on integrating we see that 

lYa.N a {z,C) = —9 (+ local constant), 

i.e., that the analytic completion of the Neumann function behaves essentially as a 
symmetric 0-connection. 

Unfortunately, the choice a = k is allowed only for simply connected domains 
because in the multiply connected case it violates (|9.1[) . In fact, 



/ Kds = 2tt(1 - g). 

JdQ. 
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On the other hand, we are allowed to take a = for any g / 1 and then 
(1 — g)T a (z, Qdz becomes a symmetric meromorphic affine connection on Q: 

(9.17) Im ((1 - g)r„(«, Qdz) = -Kds 

along dft. Moreover, it is an affine connection with simplest possible pole structure: 
simple poles with residue 1 — g at each of the points z = £ and z = £. For the 
Neumann function we get 

(9.18) Im (1 - g)W Q (z, C) = -6 (+ local constant) 
on dil. 

For the case g = 1 the situation is actually better, because in this case there 
exists a holomorphic affine connection (in agreement with the fact that in the genus 
one case the Riemann surface admits an affine structure). Just let v be the regular 
harmonic function in with Neumann boundary data ^ = — ac. This is consistent 
since fg^nds — when g = 1, and v is determined up to an additive constant. 
Now, r(z) = 2lp is the required holomorphic affine connection. 

A slightly different point of view on N a (z, £) and T a (z, £) is taken in [5T], where 
exterior domains are considered, i.e., domains f2 C P with oo G tt. In this case 



/ K ds = -2tt(1 + g) 
Jan 



and hence it is possible, for any g, to choose 

K 



a 



1 + g 

in the definition of the Neumann function. In place of (|9.17J) , (|9.18[) one then gets 

(9.19) Im((l+g)r o O>;,C)<fe) = K<fe, 

(9.20) Im(l + g)AA a (z,C) = 9 (+ constant). 

Thus, —(1 + g)T a (z, Qdz is now a symmetric affine connection on f2. (In [51j 
this connection is denoted T(z, Qdz.) However, it has not simplest possible pole 
structure. Besides the pole with residue 1 + g at z = £ it has on Q also a pole with 
residue — 2 at z — oo. Similarly on the backside fl. 

To discover the pole at infinity one has to introduce a regular (holomorphic) 
coordinate there, for example w = ~. Then 

d dw dz dw 

{w, z\\dz = — log ——dz = —2 — = 2 — . 

dz dz z w 

By (|7.4p the representative of — (1 + g)T a (z, Qdz with respect to w is 

-(i + g )r a (i,cM-)-2— , 

www 

which has residue — 2 at w — (the first term is regular at w = 0). 

As an example, consider the exterior disk: = {z 6 C : \z\ > 1} U {oo}. For 

this we have, with a = — rr- = 1- 

' i+g 

N a (z, = - log \z-C\~ log |1 - z(\ + 21og \z\ + 2 log |fl, 
(1 < \z\, |C| < oo) which gives, for the affine connection —(1 + g)T a (z, Qdz, 

„ , „ , dz dz dz 

-T a (z, Qdz = + r - 2—. 

z — Q z — i z 

In terms of w = - the same connection is given by 

_ , 1 1 . „dw dw dw n dw 

-r a -,CM- - 2— r + 7 - 2—, 

w w w w — £ w — C w 
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valid actually in all fl (\w\ < 1). Here one sees clearly the pole of residue —2 at 
w = (i.e., z = oo), along with the pole with residue one at w — ^ (i.e., z = £). 

10. Weighted Bergman spaces 

We have discussed so far the Dirichlet space A(fl) and the Bergman spaces B(il) 
and B e (n), and remarked that A(£l) and B e (fl) are connected via the simplest Bol 
operator, namely Ai : F i— > dF. In this section we shall extend the discussion to 
a sequence of weighted Dirichlet spaces Am(Q) and Bergman spaces B m (£l) and 

Be,m (^) ■ 

Let 

LU(z) 

be the Poincare metric of f2, characterized by its constant curvature and being 
complete: 

^Gauss 4, 

uj = on dfl. 

As discussed previously it gives rise to a holomorphic projective connection q(z) 
and a projective structure, namely that for which the projective coordinates are the 
coordinate on the universal covering surface of f2 when this is taken to be a disk 
or a half-plane. Let A m be the corresponding Bol operator, given in a projective 
coordinate t by A m = J^r- 

For a € C with Re a > 0, and /, g analytic in fl, smooth up to dQ, define 

(10.1) (f,g) a = f fgco^dxdy. 

Jq 

When a > is real this is an inner product on a space of analytic functions, a 
weighted Bergman space (taken to be complete, hence a Hilbert space). We denote 
this space B a (£l). It will mainly be used for a = m an integer. Clearly, Bi(il) 
is the ordinary Bergman space, and it is well-known that as a — > one gets the 
Hardy space, with the Szego inner product given by 

lim a(f, g) a = - [ fg\dz\ 
n Jan 

Remark 10.1. Our notations for weighted Bergman spaces deviate from some stan- 
dard notations used in the case of the unit disk. In for example [37] the notation A 2 a 
is used for the Bergman space with weight ^^-ui a dxdy, in place of our ui a ~ 1 dxdy, 
in the definition (jlO.ip of the inner product. In the context of automorphic forms 
the inner product (f,g) a is known as the Petersson inner product. 

If F, G are analytic in a neighborhood of Q and are kept fixed, then the inner 
product, regarded as a function of a, 

a^{F,G) a , 

is analytic for Re a > and has a meromorphic extension to all C, with simple poles 
at a = 0, —1, —2, .... The proof of this is implicit in the proof of Proposition llO.il 
below. We shall then extend the inner product (F, G) a to such values of a by 
setting 

(F, G)- m = Res Q= _ m f FGu a - l dxdy, 
Jn 

777. = 0,1,2,.... For example, for m = 0we retrieve the Szego inner product. 

Now, the main issue is that, for any m = 0, 1, 2, . . . , (F, G)— m is an inner product 
on a generalized Dirichlet space A m (fl) of analytic functions and that the Bol 
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operator is an isometric isomorphism of A m (fi) onto the subspace B e m (Q) of 'exact' 
differentials in B m (Q.)): 

A m : A m (tt) -» Be,fn(n). 

In order to elaborate the details of the above we need certain (in principle known) 
integral formulas for the A m . 

Proposition 10.1. With F, G (regarded as ^—^-forms) and g (regarded as a 
1+ 2 m -form) analytic in Q, smooth up to d£l, we have 

(10.2) f (A m F)gw m - 1 dzdz = i m - 1 (m-l)\ [ Fg(dz) 1 ^ (dz)^ , 
Jn Jon 

r _ C_ ;)m+i r 1 _ m 1+ 

(10.3) Res Q= _ m / FGu a - 1 dzdz = ^- L - l — / Fk m G{dz) — {dz)~ . 

Jn ml J an 

In other words, 

(10.4) (A m F,g) m = — -/ Fg(dz)— {dz)~ , 

z ion 

(10.5) (F,G)_ TO = t^ / FA^G^z) 1 ^^) 1 ^. 

2m! Van 

Equation \10.2)) is valid only for m — 1,2, ... , while M0.3\) , {10.$ , M0.5\) hold also 
for m = (wii/i 0! = 1). 

By choosing g = A m G in (|10.4p we obtain 

Corollary 10.2. For m = 0,1,2, ... , 

(10.6) (F, G)- m = (A m F,A m G) m , 
hence (F, F)— m > with equality if and only if A m F = 0. 

Proof. The proof of the proposition is based on ideas of Jaak Peetre, and parts of 
it have previously been outlined in |24j . 

Let {(fj} be a smooth partition of unity (i.e., tpj = 1) on such that each 
individual ipj has support within the domain of definition of a projective coordinate 
tj . In such a coordinate w is of the form Uj = a\tj | 2 + + btj + c, with a, c real and 
|fe| 2 = ac+1 (to have curvature k = —4). After an additional Mobius transformation 
we may even assume that Uj = 2 Im tj and that tj takes values in the closed upper 
half plane. Note then that dfl necessarily is mapped into the real line (since u> = 
on <9f2). 

In the coordinates tj, the coefficients F, G, g correspond to (say) Fj, Gj, gj (i.e., 

1 m 1 — m f$m 1 

F{dz)~ ~ — Fj(dtj)^~ etc.) and A m becomes On setting tfjj = tpj o tj we 
have 

f A m Fgu m - x dzdz = V f A m Fyjgu m - x dzdz 
Jn ■ Jn 

= E / -^-fiijg-jpimtjr-idtjdtj. 

j JlmtjX) m ] 

After an m-fold partial integration this becomes 

E^ 1 )" 1 / Fj^Wjg-j&lmtjr-^dtjdtj 

j k=0 3 
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Here all the boundary terms except the last one (with k — m — 1), and with 
J7F = g t m-i acting only on (2Imtj) m ~ 1 , vanish because lint,- = on R. Moving 

3 j 

back to the z-plane and taking covariance (Lemma 17. ip into account, the above 
expression becomes 

V(-l) m I FA m (i Pj goj m - 1 )dzdz + i m -\m-iy.Y' [ F^g(dz) 1 ^ 1 (dz)^ 21 

3 n 3 9Q 

= (-l) m f FgA m (u m - 1 )dzdz + i m - 1 (m-l)\ [ Fg(dz) 1 ^- (dz) 1 ^ 
Jn Jan 



l (m-l)! / Fg{dz) — {dzY 
Jau 

This proves (|10.2p . A different proof is given in [53]. 

Now we turn to (|10.3[) . In order to show that J a FGu> a ~ 1 dzdz has a meromorphic 
extension and to compute Res a =-m Jq FGuj a ~ 1 dzdz we first observe that, with 
partition of unity and notations as above, 

171 ftk sym—k 

d (E(- 1 )^(^( 2Im ^) Q+m )^rr(^ G " J )^) 

fe=0 3 
o?n-|-l £)?n-f-l 

= F j -— TT (^G j )(2Imt j r +m dt j dt j + (-i r ^-^(F J (2lmt j r+ m )i Jj G J dt j dt j . 
dtj dtj 

Therefore, 

fim—k 



'"■ r flk f\ra~k 

= / F ] -— 1 {^G ] )(2l™t j r +m dijdtj 



-(-i) m+1 (a +m) •• • (a + l)a / F J G J ip j {2lmt j ) a - 1 dt j dt j . 
The boundary terms vanish whenever Re a > and it follows that for such a, 

/ F j G j ip j (2Imt j ) a - 1 di j dt j 

Jim tj >0 
■rn+1 r »+l 

= 7 ; ? r- / F.j — n (ib i Gi)(2Imt i ) a+n dt i dt i . 

(a + m)-(« + l)ai Imt3> „ 3 c^ +lW 3>K J > 3 3 

Here the integral in the right member is an analytic function in a for Re a > — m— 1, 
hence the left member has a meromorphic extension to this range, with the residue 
at a — — m given by 

\a-lj7 JJ. _ V V / 

<9<!" +1 



Res a= _ m / F j G j il> j {2imt j ) a - , -ab j db s = ~ ' / F-^^^G^dtjA 



3 



Putting the pieces together it follows that also Res a= _ m J Q FGui a 1 dzdz has a 
meromorphic continuation to Re a > — m — 1 with residue at a = -m given by 



Res Q= _ m / FGu a 1 dzdz 



^Res Q= _ m / FjG j ip j (21mt j ) a - 1 dtjdt j 

j Jlmtj>0 
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(—i\' m + 1 r P) m 

m! j JlmtjX) ot j 
_ (-i) m+1 ^ A <9 r ' 



m! 

(_i)m+l r » 



V/ FA m { Vj G)(dz) — {dz)- 
' J an 



to! 

3 



(-i) m + 1 , 

^ " ' FA m G(dz) — (dzy 



This finishes the proof of the proposition. □ 



The corollary shows that the hermitean form [F,G)- m is positive semidefinite. 
Set 

A m (SY) = {F analytic in ft : (F, F)_ m < oo}. 
The spaces B m {^L) : with inner product (f,g) m , were defined after (|10.ip , and by 
(|10.6p the Bol operator A m is an isometry 

A m : A m (Q) -> B m {Q). 

However, this map is neither injective nor surjective when to > 0. We shall say 
something both about its kernel and its cokernel. 

As to the kernel, recall that A m in terms of any projective coordinate t is simply 

qtix 1 m 1 

A m = With F(z)(dz)^~ a holomorphic differential of order ^ LL 1 it follows 
that A m F = if and only if F becomes a polynomial of degree < to — 1 when 
expressed in any projective coordinate. Accordingly, we denote by 

P m (ty = {F e A m (n) : F holomorphic in n, A m F = 0} 

the kernel of A m . The image of A m is by definition 

£ m , e (fi) = {A m F G B m (Q) : F G A m (fl)}. 

Then we have the exact sequence 

-> P m (0) -> A m {Q) ^ B m , e {Q) -> 0. 

In other words, 

and (•, -)_ TO is positive definite on the quotient space A m (Q)/P m (£l). 

Instead of having (-,-)_ TO defined on a quotient space of A m (Q) it might be 
desirable to have it defined on a subspace of A m (Vl). One advantage then is that 
the elements in the space become functions, hence it will be possible to discuss 
reproducing kernels. There seems to be no canonical choice of such a subspace, but 
in principle it is obtained by imposing normalization conditions. To this purpose, 
choose a linear operator 

a : A m (fi) -► C m 

which does not degenerate on P m (f2), i.e., such that a(P m (f2)) = C m . Then we 
define 

A m (Cl) = {FeA m (Cl): a {F)=0} 
(the dependence on a is suppressed in the notation). Then A m (Q) = A m (p,)/P m (Q) 
and A m is an isometric isomorphism 

Am . ^4 m (ll) > _s e m (r^), 

as desired. 
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In the special case m = 1 we may choose a on the form 

a(F) = / Fads 
Jen 

for some boundary function a with J gn ads ^ and then we retrieve the previously 
discussed (see after (19. 14p ) Dirichlet space, i.e., Ai(fl) — A(fl). 

The cokernel B m {VL) / B m<e (p) of A m can be identified with the orthogonal com- 
plement of B mtf ,(yi) in B m (Q), which we write simply as B m ^ e {0)^ . And for this 
finite dimensional space (see [24] for the dimension) we have the following descrip- 
tion. 

Proposition 10.3. Bm.^Vl) 1 - consists of those elements in B m (U) which extend 
to the Schottky double O as holomorphic differentials of order ™ ■ 

Proof. By definition, g £ B mie (f2)- L if and only if 

(A m F,g) m = for all F G A m (Q), 
Using (|10.4p this becomes 

/ Fg(dz) k ^{dz)^ rL =0 
Jdn 

for all F £ A m (yi). In terms of the function 

(10.7) h(z) = ~gJz)(dz)= 1 ? !L (dz) 1 ^ =-gJz)T(z)- 1 - m , 

defined on dil, this becomes 

/ Fhdz = 
Jan 

for all F S A m (p,), which implies that h has a holomorphic extension to f2. This 
proves the proposition because h represents the continuation of g to the back-side 
of the Schottky double. Note that (|10.7[) is an instance of f2~5i □ 

Example 10.1. The case m = 1 is the well-known [55] fact that Bi, e (0) x consists 
of the abelian differentials of the first kind. 

11. Reproducing kernels 

11.1. Reproducing kernels for weighted Bergman spaces. Each of the weighted 
Bergman spaces B m (fl) and B m ^ e (£l), m > 0, have reproducing kernels K m (z,(), 
K m ,e{z,Q^ which extend to the Schottky double as differentials of order in 
each variable. The continuations to the backside are represented by the adjoint 
kernels L m (z,C), i m , e (z,C), which have singularities 

L m (z, C) = —7 ttzttt + less singular terms 

ir(z — Q m + L 

(similarly for L mje (z,£)). See [21] for details and proofs. Recall also that Ki(z,() 
is the ordinary Bergman kernel and Kq(z, £) the Szego kernel. 

For the spaces A m (Q) and A m {ty the situation is not quite that good. But 
at least A m (^) is a Hilbert space of functions for which all point evaluations are 
continuous linear functionals, and hence it has a reproducing kernel, which we 
denote IC m (z 7 £). This should be thought of as a differential of order in each 
of z and £, but its extension to the Schottky double is cumbersome because of 
appearance of branch points and multi-valuedness on the backside, as have already 
been observed in the case of Dirichlet space, m = —1, Example 19.11 Still it is 
possible to define indirectly a kind of (multi-valued) adjoint kernel, which we denote 
£ TO (z,C). In fact, we simply define C m (zX) to be any solution of equation (|11.2[) 
below. 
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Theorem 11.1. Let]C m (z,C) denote the reproducing kernel for A m (Q) and C m (z,£) 
the adjoint kernel. Then, 

(11-1) A m A m /C m (z, () = K mte (z,(), 

(11-2) A m A m C m (z,0 = L mie (z,C), 

where the first A m acts on the z-variable and the second one on The leading 
term in the singularity of C m is given by 

jm— 1 

£ m (z, C) = T7 ttt(z - C) m_1 log(z - C) + less singular terms. 

irm\{m —1)1 

Proof. The formula (jll.ip follows immediately by letting A m act on both members, 
as functions of C, in the defining equation (reproducing property) 

F(C) = (F,/C m (-,C))-m (FeAmin)) 

for JC m (z, C) and then applying (|10.6|) . 

As (|11.2I) was simply taken as the definition of £ m (z, £) it just remains to prove 
the form of the singularity. This is a matter of computation, which we omit. One 
has to check that the claimed singularity for £ m (z, C) match with that of L m ^ e (z, Q 
(which is the same as that of L m (z,£)) under pi.2p . The computation may be 
performed in a projective coordinate because the assertion only concerns the leading 
term of the singularity, which is the same in any coordinate. □ 

Example 11.1. For the unit disk we have 

•m— 1 nt— 1 m— 1 

>C m (z, = f, 777(1 - zQ m ~ l log(l -zQ+Yl + E o*(0* fc ' 

irm\(m — 1)! ' 

y ' fc=0 fc=0 

•m— 1 771—1 m— 1 

o = — f, — 777(2 - cr- 1 io g (2 -o + E a ^ z )c fe + E a *(0* fe > 

7rm! m — 1)! ' 

v ; fe=o fc=o 

where the analytic functions au depend on the normalization chosen in the definition 
of A m (fl). If the normalization for example is that F(0) = • • • = -F,„_i(0) = for 
F € A m (£l), then all the aj are zero and the kernels consist of just the first term. 

11.2. The reproducing kernel as a resolvant. Interchanging the roles of F and 
G in (|10.5|) gives, after conjugation, 

(F,G)_ m = -t^ / (A m F)G(dz) (dz) . 
2m! Jon 

With F e .A m (fi) and G = K m (z, () we get 

(11.3) F (() = - i -^ [ A m F(z)£ m (z, C) (dz) ^ (dz) ^ . 

2m! Jan 

This formula says that K, m represents the inverse, or resolvant, of A m , in terms of 
the boundary values. In the case m = 2, (111. 3f> becomes 

(11.4) F (0 = ] [ A 2 F(z)lC^z~0(dz^(dz)-^ (Cefl), 

4 Jan 

which resembles a formula of R. J. V. Jackson [313] . 
Example 11.2. For the unit disk (| 1 1 .4[) becomes 

F(C) = tt^t / F"(z)(l - 20 log(l - zQzdz, 
2tti Jan 

if ^ 2 (B) is normalized by F{0) = F'(0) = 0, (i.e., a(F) = (F(0), F'(0))). 



CRITICAL POINTS OF GREEN'S FUNCTION AND GEOMETRIC FUNCTION THEORY 57 



Equation (| 1 1 .4[) is valid for £ £ f2. For £ ^ it does not make sense because, as is 
seen clearly in Examples 19.11 and 111.11 even though the kernel extends analytically 
across <9f2, the extended kernel is not single-valued in any neighborhood of f2. 
However, it is possible to let C G ^ approach dil from inside and get a sensible 
version of (|11.4p for ( £ dil. This is what Jackson [3D] does. The equation (|11.4[) 
then solely deals with objects defined on dfl and it is appropriate to consider the 
kernel lC 2 (z, C) as a resolvent of A2. This kernel is not smooth on the boundary, as 
is seen in the example of the unit disk, there is a discontinuity in the first derivative 
caused by the jump in the imaginary part of log(l — z£) as z passes 

Remark 11.1. Since IC 2 (z, () is to be considered as a form of degree —1/2 in each of z 
and £ it is natural to compare it with 1/Kq(z, £), 1/ y/K\(z, C) and, on the diagonal, 
compare all these with u>(z). Certainly, 1/ttKq(z, z) = 1/ \JirK\(z, z) — ui(z) in 
the simply connected case. However, £2(2, C) is of a different nature because its 
continuation £2(2, C) to the back-side of the Schottky double is not single- valued. 
Similar remarks as above apply to the equation (|11.3[) for all values of m. 

12. The prepotential of a second order linear DE 

12.1. The method of Faraggi and Matone. Here we shall discuss a further 
topic related to projective structures. Consider a differential equation of the form 
(17. lip in general. We write it as 

(12.1) u" + ]^Q{z)u = 0, 

where Q(z) is a holomorphic function, say in a neighborhood of z = 0. Note that 
the Wronskian 

(12.2) W = W(ui, u 2 ) = uw'2 - u[u 2 

of any pair u\, u 2 of solutions is constant. If ui(0) 7^ 0, then any other solution u 2 
is obtained from u\ by 

(12.3) u 2 (z) = Ul (z)( f Z W£L + C) 



for suitable constants W (which then becomes the Wronskian) and C. 

An interesting approach to the problem of producing further solutions from a 
given one has been considered for instance by A. E. Faraggi and M. Matone [IF] . 
The authors introduce a "prepotential" F{u) = !F Ul (u), a function of a complex 
variable u, which functionally depends also on the first solution u\ (alternatively, 
depends on Q if initial conditions for u\ are specified), to the effect that a second 
solution U2 is obtained by taking the derivative with respect to u: 

As noted in [16] . !F(u) is essentially a Legendre transform, namely of the indepen- 
dent variable z considered as a function of any projective coordinate t associated to 
the projective structure given by Q(z). We proceed to explain briefly these issues, 
going slightly beyond [TFJ . 

Fixing a value W 7^ of the Wronskian, consider pairs u\ = U\{z), u 2 — u 2 (z) 
of solutions, holomorphic in a neighborhood of z — and subject to (112. 2 [I . i.e., 

(12.5) u\du 2 — u 2 du\ — Wdz. 
We assume that ui(0) 7^ and set 

(12.6) s = u\, t=—. 

Ml 
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Then 

~r = — , {t,z} 2 = Q{z), 

dz s 

and also, for example, {z,t}i — • In terms of s and i, (|12.5p becomes 

(12.7) sdt = Wdz. 

We note also that, in terms of s and t, the (m— l)-fold symmetric product S' m_1 (L) 
of the operator L = -j^ + \Q{z), mentioned briefly after Lemma [7.11 will have 

solutions generated by s^~t k , k = 0, 1, . . . , m — 1. 

In a neighborhood of z = we can invert i = t(z) to consider 2 as a function of 
t: z — z(t). Then =| = Assuming for a moment that ^ we can form the 
Legendre transform of Wz(t). It is 

(12.8) C(s) = st- Wz(t) with t chosen so that s = W—Q. 

dt 

Note that the final equation assures that the variable t keeps its meaning as t — ^ 
when s = u\. 

By a direct computation, or by using that the Legendre transform is involutive, 
one realizes that 

dC(s) _^ 
ds 

with s, t related to u\, u 2 as above. This shows that 

1 ds dC(s) IdLiu?) 
u 2 =u 1 -t= — = ^-ii, 

2 du\ ds 2 dui 

hence establishes the Legendre transform as essentially the desired prepotential: 

Hu) = \c{u 2 ). 

An alternative and slightly more general approach, which makes sense also if 

,2 

^# = 0, is to consider u\, u 2 and z as independent variables, or coordinates, in 
a three dimensional space. Then (|12.6p is simply a coordinate transformation (in 
two of the variables) and (| 1 2 . 5[) , (|12.7|) should be thought of as defining a contact 
structure (see, e.g., [2J, Appendix 4). With 

(12.9) C = st-Wz = Ul u 2 - Wz, 

now considered as a function in the three dimensional space, we have 

dC = tds = 2u 2 du\ 
when the contact structure is taken into account. This gives again 
_ dC 1 _ 1 dC I 

\sdt—Wdz^ 2 dll\ ^ U l^ u 2—U2dui—Wdz' 

To derive explicit formulas for T(u) it is convenient to assume that the first 
solution i*i is invertible near z = 0, i.e., to assume that u[(0) ^ 0. Since the 
independent variable u in T(u) will finally, in (|12.4[) . be assigned to have the value 
ui(z), it is natural to invert this relation to z = (u). Substituting into (|12.3[) . 
(|12.9p and recalling that T is half of £ one arrives at 

1 1 /• u r 1 (") wdc 

(12.10) T U1 (u) = -( Ul u 2 - Wz) = -(u 2 + Cu 2 - Wui\u)). 

After a partial integration and a change of variable this gives the formula presented 
in [IS]: 

(12.11) = u 2 ( r + %) = Wu 2 f U ^ du 
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Here the final integral is an "indefinite integral". One readily verifies that (|12.4[) 
indeed holds. 

As a final issue one may notice (as in [16p that the prepotential satisfies a third 
order differential equation. This is obtained by considering u 2 as a function of u\ 
(in place of z). This renders (|12.ip of the form 

d 2 u 2 Q , du 2 , 3 

Replacing u\ by u and u 2 by dJ ^^ gives 

d a F(u) _ _Q_ d 2 T{u) dT(u) 3 
du 3 ~ ZW 2 ^ du 2 cfa~ ' ' 

12.2. Example. Consider the differential equation 

u" + u = 0, 

i.e., (|12.1[> with Q = 2, and choose 

iti (z) = cosz. 

Then 

f z Wd( 

u 2 (z) — cos z( 7T- + C) = W sin z + C cosz. 

Jo cos C 

Using l|12.10p the prepotential becomes 

^ /-arccos it ^J/^ 

T u Au) = - (u 2 1 5— +Cu 2 —W arccos u) = — (uyl — u 2 — arccos u) +— u 

2 Jo cos^ C 2 2 

Somewhat surprisingly perhaps, we start with a very simple differential equation 
and arrive at a rather complicated prepotential. Note that u\{z) = cosz violates 
the assumption of being invertible at z — 0. This causes J- Ul (u) to have a branch 
point at u = Ui(0) = 1. 

13. Glossary of notations 

• B(a,r) = {z G C : |z - a\ < r}, D = D(0, 1). 

• P = C U {00}, the Riemann sphere. 

• A rifi = {z G C : r < |z| < R}. 

• usually denotes a finitely connected domain in P, with boundary com- 
ponents denoted ro,...r g , g > (each Tj consisting of more than one 
point). 

• Q = O U d£l U f2, the Schottky double of O, a symmetric compact Riemann 
surface of genus g. 

• J : Q — > fi, the anticonformal involution on fi. 

• ay, /3j (J = 1, ...,g), canonical homology basis on a compact Riemann 
surface of genus g. 

• d(z, A) = dist (z, A), distance from a point z to a set A; d(z) = d(z, dil) if 
2 G SI 

• *a->, the Hodge star of a differential form uj, for example *{adx + bdy) = 
—bdx + ady, *dz — —idz. If u is a harmonic function then *du = d(u*) 
where u* is a harmonic conjugate of u. 

• V(z, W, a,b), a fundamental potential on a compact Riemann surface. 

• E(z,(), the Schottky-Klein prime function. 

• G(z, C), the Green's function of a domain fi C P. 

• H(z, C), the regular part of the Green's function 

• G(z, () = G(z, ()+iG*(z, (), the analytic completion of G(z, () with respect 
to z. Similarly for TL{z,Q). 
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• P{z,Q) — Qn ' * ne Poisson kernel outward normal derivative). 

• Gj(z, (), the hydrodynamic Green function with circulations 7. 

• Gj(z, £), the analytic completion of G 7 (z, £) with respect to z. 

• N a (z, £), the Neumann function with Neumann data —a on dtt. 

• N a {z, C), the analytic completion of N a (z, () with respect to z 

• v a -b, the abclian differential of the third kind with poles at a, b and having 
purely imaginary periods. 

• Ld a -b, the abelian differential of the third kind with poles at a, b and with 
vanishing Oj-periods. 

• u) a -b, the abelian differential of the third kind with poles at a, b and with 
vanishing /^-periods. 

• A m (£l), A(fl) — Ai(£l): spaces of analytic functions in £1 (to = 0, f , 2, . . . ). 
Essentially weighted Bergman spaces of negative index. A (Q) is Hardy 
space. 

• A m (ty, A(Cl) — „4i(£l): subspace of A m (n) defined by a normalization (to 
make the inner product positive definite). A(£l) is Dirichlet space. 

• B m (il), B(il) = -Bi(f2): Weighted Bergman spaces of positive index to = 
1,2, ... . B(fl) ordinary Bergman space. 

• B e (Q), B m ^ e (Q) (to = 1,2, . . .), the subspaces of B(Q), B m (il) consisting 
of "exact" differentials with respect to A m . 

• P m (0) = {F e A m (n) ■ A m F = 0} = {Fe A m {Q) : (/,/)_ m - 0}, the 
space of — j^:s order differentials in f2 which expressed in any projective 
coordinate are polynomials of degree < to — 1 (to =1,2,...). 

• H(fl), H e (£l): spaces of harmonic functions. 

• {fi9)mi the inner product on B m (fl) for m > 1, on A- m (ty for m < 0. 

• D(f,g), Dirichlet inner product. 

• IC m (z, £), reproducing kernel for A m (£l); tC(z, £) = K\{z, (). 

• £ m (z, C), adjoint kernel for A m (Sl); £(z, () = d(z, £). 

• K(z, () = — ^ d Q^df* ' ^ e Bergman kernel, reproducing kernel for B(fl). 

• K e (z,() = — ^ 9 ^dc'^ ' ^ e re( i uce d Bergman kernel, reproducing kernel 
for B e (Q). 

• K m (z, £), the reproducing kernel for B m (Vl). 

• K m , e (z,(), the reproducing kernel for B m , e (Q). 

• L(z,Q = f d Qzd'c^ ' ^ e Schiffer kernel, or adjoint Bergman kernel. (Simi- 
larly for L m (z, C), L m ^ e (z, Q.) 

• £(z, C), the regular part of L(z, (). 

• k(z,(), k e (z,Q, reproducing kernels for H($l) and H e (Q) respectively. 

• M(z, £), the ordinary Martin kernel. 

• F(z, (), the Martin kernel for the gradient structure. 

• ds, arc-length diffcrental along a curve. 

• k, the curvature of a curve in the complex plane (and sometimes short for 

^Gauss) ■ 

• da = p(z)\dz\ — ^J0y = e p ( z )\dz\, a hermitean metric, e.g., the Poincare 
metric. 

• KGauss, the Gaussian curvature of a hermitean metric. 

• S(z), the Schwarz function of an analytic curve (S(z) = z on the curve, 
S(z) analytic in a full neighborhood of the curve). 

• T{z) = the unit tangent vector on a curve (and its analytic extension 
to a neighborhood of the curve, if the curve is analytic). 

• {z, t}k, a differential expression appearing in the definition of a fc-connection 
(fc = 0, f , 2). (The Schwarzian derivative if k = 2.) 
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• Vfe, the covariant derivative for an affine connection, when acting on fc:th 
order differentials. 

• A m , the m:th order Bol operator, the covariant derivative for a projective 
connection, acting on differentials of order 

• V, the gradient. 

• C, the Legendre transform. 

• 5"™, the m-fold symmetric product (of a differential operator). 



References 

[1] L. V. Ahlfors: Conformal Invariants, Topics in Geometric Functions Theory, McGraw-Hill, 
New York, 1973. 

[2] V. Arnold: Mathematical Methods in Classical Mechanics, Springer- Verlag, New York, 1989. 
[3] S. Bell: The Cauchy Transform, Potential Theory, and Conformal Mapping, Studies in 

Advanced Mathematics, CRC Press, 1992. 
[4] S. Bergman: The Kernel Function and Conformal Mapping, American Mathematical Society, 

1950. 

[5] L. Bieberbach: Au = e u und die Automorphen Funktionen, Math. Ann. 77 (1916), 173-212. 
[6] G. Bol: Invarianten linearcr Diflcrcntialglcichungen, Abh. Math. Sem. Univ. Hamburger 
Univ., 16(1949), 1-28. 

[7] C. Constantincscu, A. Cornea: Ideale Rdnder Riemannscher Fldchen, Springer Verlag, 
Berlin, 1963. 

[8] R. Courant, D. Hilbert: Methods of Mathematical Physics, Intcrscicnce Publishers, New 
York, volume 1, 1953. 

[9] D. G. Crowdy, J. S. Marshall: Green's functions for Laplace's equation in multiply connected 

domains, IMA J. Appl. Math., 72 (2007), 1-24. 
[10] P. Davis, H. Pollak: A theorem for kernel functions, Proc. Am. Math. Soc, 2 (1951), 686-690. 
[11] P. Davis: The Scharz function and its Applications, Carus Mathematical Monographs, 1974. 
[12] B. Dubrovin: Geometry of 2D topological field theory, Integrable systems and quantum 

groups, M. Francaviglia and S. Greco eds., Springer Lect. Notes in Math. 1260, 1996, pp. 

120-348. 

[13] G. Duff, D. Naylor: Differential Equations of Applied Mathematics Wiley & Sons, New York, 
1966. 

[14] M. Eichler, D. Zagier: On the zeros of the Weierstrass p-function, Math. Ann., 258 (1982), 
399-407. 

[15] B. Epstein: The kernel function and conformal invariants, Journal of Mathematics and Me- 
chanics, Vol 7, No. 6 (1958). 

[16] A. E. Faraggi, M. Matonc: Duality of x and vb an d a Statistical Interpretation of Space in 
Quantum Mechanics. Phys. Rev. Lett. 78 (1997), 163 . 

[17] H. M. Farkas, I. Kra: Riemann Surfaces, Springer Verlag, 1980. 

[18] J. Fay: Theta functions on Riemann surfaces, Lectures Notes in Math. 352, Springer Verlag, 
1973. 

[19] M. Flucher: Variational Problems with Concentration, Birkhauser, 1999. 
[20] T. Frankel: The Geometry of Physics. An Introduction. Second edition, Cambridge University 
Press, Cambridge, 2004. 

[21] I. S. Gal: The Martin boundary, Proceedings of the National Academy of Sciences of the 

United States of America, 48 (1962), 775-780. 
[22] R. C. Gunning: Lectures on Riemann Surfaces, Princeton University Press. Princeton, 1966. 
[23] B. Gustafsson, J. Pcctrc: Notes on projective structures on complex manifolds, Nagoya Math. 

J., 116 (1989), 63-88. 

[24] B. Gustafsson, J. Peetre: Hankel forms on multiply connected plane domains. Part two. The 

case of higher connectivity, Complex Variables, 13 (1990), 239-250. 
[25] B. Gustafsson, A. Vasilev: Conformal and Potential Analysis in Hele-Shaw Cells, Birkhauser, 

Basel, 2006. 

[26] N. S. Hawley, M. Schiffer: Half order differentials on Riemann surfaces, Acta Mathematica, 
115, 199-236, 1966. 

[27] H. Hedenmalm, B. Korenblum, K. Zhu: Theory of Bergman Spaces, Springer- Verlag, New 
York, 2000. 

[28] D. A. Hejhal: Theta Functions, Kernel Functions and Abelian Integrals, Memoirs of the 
American Mathematical Society, 1972. 



62 BJORN GUSTAFSSON AND AHMED SEBBAR 

[29] D. A. Hejhal: Some remarks on kernel kunctions and Abelian differentials, Arch. Rat. Mech. 

Anal., 52 (1973), 199-204. 
[30] R. J. V. Jackson: On the boundary values of Riemann's mapping theorem, Trans Amer. 

Math. Soc. 259 (1980), 281-297. 
[31] R. P. Jerrard, L. A. Rubel: On the curvature of the level lines of a harmonic function, Proc. 

Am. Math. Soc, 14 (1963), 29-32. 
[32] J. L. Kelley: General Topology, Van Nostrand, 1955. 

[33] N. Koblitz: Introduction to Elliptic Curves and Modular Forms, 2nd ed., Spriger- Verlag, 
New York, 1993. 

[34] I. K. Kostov, I. Krichever, M. Mineev-Weinstein, P. Wiegmann, A. Zabrodin: r-function for 

analytic curves, in Random matrices and their applications, MSRI publications 40 Cambridge 

Academic Press, Cambridge, 2001. 
[35] I. Krichever, A. Marshakov, A. Zabrodin: Inegrable structure of the Dirichlet boundary 

problem in multiply-connected domains, Comm. Math. Phys. 259 (2005), 1-44. 
[36] T. Kubo: On the potential defined in a domain, Proc. Japan Acad., 25,(1949) 123- 125. 
[37] C. C. Lin: On the Motion of Vortices in Two Dimensions, University of Toronto Press, 

Totonto, 1943. 

[38] P. Levy: Sur les valeurs de la fonction de Green dans le voisinage du contour, Bull. Sci. Math. 
34 (1910), 1-4. 

[39] A. J. Maria: Concerning the equilibrium point of Green's function for an annulus, Duke 

Math. J., 1 (1935), 491-495. 
[40] J. Milnor: Dynamics in One Complex Variable, Introductory Lectures, 2nd edition, Vieweg 

Verlag, Weisbaden, 2000. 
[41] M. Mineev, P. Wiegmann, A. Zabrodin: Integrable structure and interface dynamics, Phys. 

Rev. Lett. 84 (2000), 5106-5109. 
[42] D. Mumford: Tata Lectures on Theta II, Birkhauser, Basel, 1984 

[43] P. K. Newton: The N -Vortex Problem. Analytical Techniques, Springer- Verlag, New York, 
2001. 

[44] H. Poincare: Les fonctions fuchsiennes ct l'cquation Am = e u , J. Math. Pures Appi, (5) 4 
(1898), 137-230. 

[45] J. Polchinski: String Theory, Volume I, An Introduction to the Bosonic String, Cambridge 

University Press, Cambridge, 2001. 
[46] G. Polya, G. Szego: Problems and Theorems in Analysis I, Springer- Verlag, New- York, 1972. 
[47] M. Protter, H. Weinberger: Maximum Principles in differential Equations, Prentice-Hall, 

Englcwood Cliffs, 1967. 

[48] S. Ramanujan: On certain arithmetical functions, Trans. Cambridge Philos. Soc. 22 (1916), 
159- 184. 

[49] L. Sario, K. Oikawa: Capacity Functions, Springer Verlag, 1969. 

[50] M. Schiffer: The span of multiply connected domains, Duke Math. J. 10 (1943), 209 - 216. 
[51] M. Schiffer, N. S. Hawley: Connections and conformal mapping, Acta Math., 107 (1962), 
175-274. 

[52] M. Schiffer, D. C. Spencer: Functionals of Finite Riemann Surfaces, Princeton University 
Press, Princeton, 1954. 

[53] A. Sebbar, T. Falliero: Equilibrium points on the annulus and Eiscntein series. Proc. Amer. 

Math. Soc. 132 (2007), 313 - 328. 
[54] H. S. Shapiro: The Schwarz Function and its Generalization to Higher Dimensions John 

Wiley&Sons, New York, 1992. 
[55] C. L. Siegel: Lectures on Advanced Analytic Number Theory, Tata Institute. Tata Institute 

of Fundamental Research lectures on Mathematics, No 23. Tata Institute of Fundamental 

Research, Bombay 1965. 

[56] A. Y. Solynin: A note on Equilibrium points of Green's Functions. Proc. Am. Math. Soc. 136 
(2008), 1019-1021. 

[57] L. A. Takhtajan: Quantum field theories on an algebraic curve, Letters in Mathematical 

Physics 52 (2000), 79-91. 
[58] L. A. Takhtajan, L.-P. Teo: Quantum Liouville theory in the background field formalism I. 

Compact Riemann surfaces, Comm. Math. Phys. 268 (2006), 135—197. 
[59] A. Yamada: Positive differentials, theta functions and Hardy H 2 kernels. Proc. Amer. Math. 

Soc. (5) 127 (1999), 1399-1408. 
[60] K. Zarankicwicz: Uber ein numerisches Verfahrcn zur konformen Abbildung zweifach zusam- 

menhangender Gebiete. Zeit. f. Ang. Math. u. Mech., 14 (1934), 97-104. 



CRITICAL POINTS OF GREEN'S FUNCTION AND GEOMETRIC FUNCTION THEORY 63 



B. GUSTAFSSON, DEPARMENT OF MATHEMATICS, KTH, 100 44 STOCKHOLM, SWEDEN. 
E-mail address: gbjornaSkth.se 

A. Sebbar, Universite Bordeaux I, Institut de Mathematiques, 33405 Talence, France 
E-mail address: sebbaraSmath.u-bordeauxl.fr 



